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ABSTRACT 


A  procedure  for  the  development  of  a  simple  boundary  layer 
turbulence  model  to  account  for  different  physical  effects  Is 
described;  the  method  Is  applied  here  to  produce  models  for  both 
pressure  gradient  and  mainstream  turbulence  effects.  Asymptotic 
theory  Is  used  to  Isolate  the  leading  terms  In  an  expansion  for 
the  mean  velocity  profile  for  high  Reynolds  nutrtjers  for  both  the 
Inner  and  outer  regions  of  a  nominally  steady  two-dimensional 
boundary  layer.  The  velocity  profile  In  the  outer  layer  satis¬ 
fies  a  partial  differential  equation  containing  a  Reynolds  stress 
term  and  this  term  Is  modeled  by  a  simple  eddy  viscosity  function 
which  contains  two  parameters.  The  velocity  profile  In  the  Inner 
wall  layer  Is  modeled  using  an  analytical  expression  which  has 
been  previously  derived  by  consideration  of  the  observed  charac¬ 
teristics  of  the  time-dependent  flow  In  the  wall  layer  and  which 
contains  a  single  Independent  parameter.  For  a  self-similar 
flow,  the  outer  layer  equation  becomes  an  ordinary  differential 
equation;  this  equation  Is  solved  numerically  and  In  conjunc¬ 
tion  with  the  analytical  Inner  layer  profile,  a  composite  pro¬ 
file  spanning  the  entire  boundary  layer  Is  defined.  This  com¬ 
posite  profile  contains  three  parameters  which  may  be  adjusted 
systematically  to  obtain  a  best  fit  to  a  given  set  of  experimen¬ 
tal  data. 
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A  computer  optimization  code  Is  described  in  which  any  or 
all  of  the  three  profile  parameters  may  be  varied.  This  optimi¬ 
zation  code  may  be  used  simply  to  obtain  a  close  analytical 
representation  of  a  given  set  of  data.  The  primary  use  described 
here,  however,  is  to  develop  correlations  for  various  physical 
effects  from  the  results  of  the  optimizations.  In  particular, 
correlations  for  the  effects  of  mainstream  pressure  gradients 
and  mainstream  turbulence  for  the  profile  parameters  are  given. 

In  principle,  these  correlations  may  then  be  used  in  a  predic¬ 
tion  method. 
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1.  INTRODUCTION 

Turbulent  boundary  layers  occur  in  a  wide  variety  of 
engineering  applications  including,  for  example,  flows  over 
turbine  blades,  airfoils,  and  in  subsonic  diffusers.  In  these 
and  many  other  situations  it  is  important  to  be  able  to  develop 
the  capability  to  accurately  predict  characteristics  such  as 
skin  friction,  lift,  drag,  and  the  onset  of  boundary  layer 
separation.  At  present,  the  complex  nature  of  turbulence 
seems  to  preclude  any  predictive  analysis  based  on  first  prin¬ 
ciples.  For  a  turbulent  boundary  layer  which  is  steady  in  the 
time-mean  sense,  the  classical  approach  to  prediction  has  been 
to  attempt  to  deal  with  the  Reynolds  time-averaged  equations 
for  the  mean  velocity  components.  These  equations  contain 
unknown  functions  referred  to  as  Reynolds  stresses;  physically 
these  terms  are  related  to  long  time  averages  of  products  of 
fluctuations  of  the  velocity  components  about  their  mean  values. 
As  a  first  step,  the  development  of  any  prediction  method 
requires  a  model  to  represent  the  behavior  of  the  Reynolds 
stress  terms.  The  objective  of  this  study  is  to  develop  an 
approach  which  can  be  used  to  investigate  the  influence  of 
various  physical  effects  on  the  turbulent  boundary  layer  and 
to  incorporate  these  effects  in  a  simple  turbulence  model. 

To  this  end,  a  general  optimization  code  is  developed  in 


which  the  parameters  in  the  turbulence  model  may  be  varied 
systematically  to  obtain  the  best  representation  of  measured 
mean  profile  data.  As  one  example  of  how  this  code  may  be 
used,  the  method  is  applied  to  determine  the  influence  of 
pressure  gradients  on  the  mean  profile  and  a  correlation  for 
pressure  gradient  effects  is  developed.  As  a  second  example, 
the  method  is  applied  to  some  recent  constant  pressure  pro¬ 
file  data  with  mainstream  turbulence  and  a  correlation  is  obtained 
to  reflect  the  effects  of  mainstream  turbulence. 

An  immediate  problem  which  arises  in  the  modeling  process 
is  that  the  dynamics  of  turbulent  boundary  layers  are  not  well 
understood.  Previous  investigators  have  attempted  to  resolve 
this  problem  by  postulating  the  functional  form  for  the  Reynolds 
stress  terms.  These  functional  forms  normally  contain  a  num¬ 
ber  of  unknown  constants  which  are  selected  in  a  procedure  often 
known  as  "computer  optimization."  The  details  of  this  procedure 
vary  with  the  originator  of  the  particular  model  and  invariably 
are  not  well  documented.  However,  the  general  approach  is 
that  particular  data  sets  usually  consisting  of  measured  velo¬ 
city  profiles  at  numerous  streamwise  locations  are  "predicted" 
using  various  combinations  of  values  of  the  "constants"  asso¬ 
ciated  with  a  given  model.  Some  type  of  subjective  judgement 
is  then  made  as  to  which  set  of  values  of  the  constants  best 
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"predict"  as  many  data  sets  as  possible.  The  approach  adopted 
in  this  study  is  rather  different  and  will  now  be  discussed  in 
detail. 

An  ideal  approach  to  the  modeling  problem  would  be  to 
isolate  the  primary  features  in  the  time-dependent  turbulent 
flow  and  to  pattern  the  turbulence  model  after  motions  which 
reflect  the  true  physics  of  the  turbulent  boundary  layer.  Such 
a  model  would  adequately  replace  the  information  lost  in  the 
time  averaging  and  would  presumeably  not  contain  or  require  a 
large  number  of  adjustable  constants.  This  approach  is  par¬ 
tially  adopted  in  this  study  in  the  model  used  for  the  inner 
region  of  the  turbulent  boundary  layer. 

Since  the  objective  is  to  isolate  which  effects  are  impor¬ 
tant  and  which  are  not  for  large  Reynolds  number,  singular  per¬ 
turbation  theory  and  the  n«thod  of  matched  asymptotic  expan¬ 
sions  are  used  throughout  the  present  study.  It  is  worthwhile 
at  this  stage  to  sunmarize  some  of  the  main  results  concerning 
the  asymptotic  structure  of  the  time  mean  boundary  layer  equa¬ 
tions  in  the  limit  of  large  Reynolds  nimbers;  these  results 
have  formally  been  demonstrated  by  Fendell  (1972)  and  Mel  lor 
(1972).  The  turbulent  boundary  layer  is  a  composite  double 
layer  consisting  of  a  relatively  thick  outer  layer  having  a 
thickness  OCu*),  (where  u*  »  u^/Ug(x)  is  the  ratio  of  friction 
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to  mainstream  velocity)  and  a  thin  inner  layer  having  a  thick¬ 
ness  0(l/(Re  i/')),  where  Re  is  the  Reynolds  number.  In  the  join¬ 
ing  region  between  the  two  layers,  the  velocity  profile  must  be 
logarithmic  for  a  self  consistent  asymptotic  description.  The 
complete  mathematical  results  for  the  streamwise  momentijii 
equations  are  sumnarized  in  §2.  One  particularly  important 
result  concerns  the  inner  layer  and  is  that,  to  leading  order, 
the  convective  terms  are  negligible  in  the  time  mean  equations; 
consequently,  if  the  mean  profile  is  known  in  the  inner  layer, 
the  Reynolds  stress  may  be  calculated  and  vice  versa.  In  the 
present  study,  rather  than  a  model  for  Reynolds  stresses,  a 
model  for  the  inner  region  velocity  profile  is  used;  this  model 
is  based  upon  the  observed  coherent  structure  of  the  wall  layer 
flow  and  will  now  be  briefly  described. 

Over  the  past  decade,  it  has  been  well  documented  (see, 
for  example,  Kline  &  Runstadler,  1959;  Kline  et  al.,  1967; 

Corino  &  Brodkey,  1969;  Willmarth,  1975)  that  there  is  a  con¬ 
siderable  degree  of  ordered  structure  in  the  time  dependent 
flow  in  the  wall  layer  of  a  turbulent  boundary  layer.  In  par¬ 
ticular,  it  is  well  known  that  there  are  two  important  phases 
associated  with  an  observed  cyclic  behavior  of  the  wall  layer 
flow.  In  the  first  phase,  if  attention  is  focused  on  a  fixed 
small  area  of  the  plate,  the  wall  layer  will  be  observed  to 
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be  in  the  quiescent  state  (Kline  et  al.,  1967)  for  a  majority 


of  the  total  observation  time.  During  this  quiescent  period, 
the  wall  layer  streaks  are  observed  with  what  appear  to  be 
longitudinal  counter-rotating  vortices  between  the  streaks,  the 
wall  layer  flow  is  relatively  well  ordered  and  no  important 
interactions  occur  between  the  wall  layer  and  the  flow  in  the 
outer  portion  of  the  boundary  layer;  in  this  state,  the  wall 
layer  may  be  regarded  as  passive.  Eventually  the  second  phase 
occurs  which  is  generally  kn(»vn  as  the  bursting  phenomenon  and 
which  is  characterized  by  a  rapid  and  violent  ejection  from  the 
wall  layer  into  the  outer  layer.  The  ejection  is  of  relatively 
short  duration  and  is  followed  by  an  inrush  of  fluid  from  the 
outer  layer;  the  streak  structure  appears  again  very  rapidly  and 
another  quiescent  period  begins.  Although  many  questions  exist 
as  to  the  causes  and  effects  of  these  and  other  subsidiary 
events,  the  gross  features  of  the  cyclic  behavior  described 
above  arenowwell  established. 

To  incorporate  such  information  in  a  prediction  method 
for  the  time-mean  flow,  it  is  necessary  in  principle  to  analyze 
a  typical  event  in  the  turbulence  and  then  assess  the  contribu¬ 
tion  of  this  event  to  the  time-mean  quantities.  To  this  end. 
Walker  and  Abbott  (1977)  argue,  by  consideration  of  the  observed 
length  and  time  scales  in  the  wall  layer,  that  during  the 
quiescent  period  the  equations  for  all  three  velocity  components 
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must  be  linear  and  of  the  heat  conduction  type  in  the  limit  of 
large  Reynolds  number;  in  other  words,  the  convective  terms 
in  the  Navler-Stokes  equations  are  negligible  to  leading  order 
whenever  the  wall  layer  flow  is  in  the  quiescent  state.  Walker 
and  Scharnhorst  (1977)  then  go  on  to  consider  all  possible 
similarity  solutions  of  these  equations  which  are  compatible 
with  theory  and  experiment;  the  similarity  solutions  correspond 
to  the  organized  motion  between  streaks  during  the  quiescent 
period  and  are  radically  different  from  the  oscillatory  Stokes 
type  solutions  which  form  the  basis  of  the  Van  Driest  (1956) 
model.  Walker  and  Scharnhorst  (1977)  compute  a  time-average 
of  the  similarity  solutions  over  a  quiescent  period  and  assess 
which  solutions  produce  an  important  contribution  to  the  time- 
mean  profile  and  which  do  not.  The  final  result  is  an  analy¬ 
tical  model  which  will  be  given  in  iZ  for  the  inner  region  pro¬ 
file  and  which  contains  a  single  parameter  S  that  is  related  to 
the  mean  period  between  bursts.  The  contribution  to  the  mean 
profile  during  the  bursting  process  and  breakdown  of  the  wall 
layer  flow  is  neglected  on  the  grounds  that  the  breakdown  is 
of  short  duration  relative  to  the  quiescent  period;  note  that 
there  are  various  theoretical  reasons  as  well  as  a  body  of 
experimental  evidence  that  verify  that  the  period  of  breakdown 
must  be  small  with  respect  to  the  quiescent  period.  On  the 
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other  hand,  the  vertical  velocity  (Walker  &  Abbott,  1977  and 
Walker  &  Scharnhorst,  1977)  is  so  small  during  the  quiescent 
period  that  there  can  be  no  contribution  to  the  Reynolds  stress 
to  leading  order  during  this  period  of  time  and  the  major  con¬ 
tribution  to  the  Reynolds  stress  must  be  made  during  the  burst¬ 
ing  process.  Consequently,  in  this  theoretical  description  of 
the  wall  layer  flow,  the  dominant  contribution  to  the  mean  pro¬ 
file  occurs  during  the  quiescent  period,  while  the  major  con¬ 
tribution  to  the  Reynolds  stress  occurs  during  the  breakdown 
phase;  both  mean  quantities  are  directly  related  to  each  through 
the  leading  order  time-mean  equations  in  the  wall  layer  (since 
the  convective  terms  are  negligible  to  leading  order). 

In  principle,  it  is  desirable  to  develop  a  model  for  the 
outer  layer  which  is  also  based  on  the  observed  dynamics  of 
the  time  dependent  flow  in  the  outer  region.  However,  the  outer 
region  problem  is  more  complex  ’*  it  appears  necessary  to  model 
the  Reynolds  stress  terms  directly  and  unfortunately  the  dynamics 
of  the  outer  layer  are  not  well  understood.  In  the  outer  region, 
large  scale  motions  are  observed  which  appear  to  be  recirculating 
agglomerations  of  numerous  smaller  scale  structures;  these 
smaller  scale  structures  have  dimensions  on  the  general  magni¬ 
tude  of  100  wall  layer  units  (y*)  and  appear  at  times  to  be 
intensely  vorticular  in  nature.  At  any  stage,  these  small 
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vorticular  structures  are  in  various  stages  of  coalescence  and 
decay  and  when  they  pass  close  to  the  wall  layer,  are  observed 
to  Induce  eruptions  of  fluid.  A  number  of  authors  have  sug- 
gessted  that  it  is  the  vortex  motion  in  the  outer  layer  which 
induces  the  wall  layer  bursting  (Nychas  et  al.,  1973;  Doligalski 
&  Walker,  1978;  Walker,  1978)  and  leads  to  the  creation  of  new 
vortex  structures  in  the  outer  layer;  a  regenerative  mechanism 
for  the  production  of  new  turbulence  has  recently  been  proposed 
by  Doligalski,  Smith  &  Walker  (1981)  and  Doligalski  (1981)  on 
the  basis  of  the  observed  unsteady  effect  of  vortices  on  wall 
boundary  layers.  These  studies  demonstrate  that  for  the  vortex 
motions  considered  (twc-dimensional  vortices  convected  in  a  uni¬ 
form  flow  and  in  e  shear  flow  and  ring  vortices  above  a  plane 
wall)  that  an  eruptia*»  of  the  boundary  layer  flow  near  the  wall 
will  occur.  However,  these  studies  are  as  yet  in  an  explora¬ 
tory  stage  and,  while  a  physical  mechanism  for  the  bursting 
is  indicated,  it  is  not  yet  possible  to  develop  a  constitutive 
model  for  the  outer  layer  on  this  basis. 

The  bursting  phenomena  is  a  complex  time-dependent  viscous- 
inviscid  interaction  between  the  inner  and  outer  layers  which 
leads  to  the  introduction  of  new  vorticity  into  the  outer  layer 
and  which  occurs  intermittently.  However,  during  these  brief 
periods  of  localized  breakdown  of  the  two  layer  structure,  the 
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majority  of  the  contribution  to  the  outer  layer  (and  inner 
layer)  Reynolds  stress  occurs.  For  this  reason,  the  simplifi¬ 
cations  which  were  possible  in  obtaining  the  wall  layer  model 
are  not  possible  and  it  appears  necessary  to  consider  a  typical 
burst  and  time  average  the  results  in  order  to  model  the  u'v* 
term.  As  discussed  by  Doligalski  &  Walker  (1978)  and  Doligalski, 
Smith  &  Walker  (1981)  the  inviscid-viscous  interaction  is  a 
formidable  theoretical  problem  at  present;  moreover,  many  other 
aspects  of  the  outer  layer  are  still  not  well  understood  and, 
at  present,  development  of  an  outer  layer  model  based  on  the 
characteristics  of  time  dependent  flow  in  the  outer  layer  is 
not  feasible. 

For  these  reasons,  a  conventional  type  of  eddy  viscosity 
model  is  used  in  the  present  study  for  closure  in  the  outer 
layer;  this  type  of  model  is  used  here  because  of  its  simpli¬ 
city  and  the  good  degree  of  success  it  has  had  in  other  pre¬ 
diction  methods.  The  eddy  viscosity  hypothesis  is  coitmonly 
termed  a  first-order  closure  scheme;  it  assumes  that  there  is 
a  functional  relation  between  the  Reynolds  stress  and  the 
mean  profile  and  further  that  the  Reynolds  stress  may  be  written 
as  an  eddy  viscosity  function  times  the  mean  velocity  gradient. 
Here  a  model  for  the  eddy  viscosity  which  is  similar  to  the 
Cebeci-Smith  (1974)  and  Mellor  &  Herring  models  (1968)  will 
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be  adopted;  this  model  Is  a  monotonically  increasing  function 
of  the  outer  variable  behaving  linearly  with  slope  k  (the  von 
Karman  "constant")  near  the  logarithmic  zone  and  approaching 
a  value  K  at  the  boundary  layer  edge.  Note  that  the  wall 
layer  profile  model  contains  k  in  addition  to  the  cycle  time 
parameter  S.  Consequently  the  present  turbulence  model  contains 
the  three  parameters  (s,k,K). 

It  is  common  practice  in  many  turbulence  models  to  assume 
constant  universal  values  for  the  parameters  appearing  in  the 
model.  Scharnhorst  (1978),  for  example,  has  attempted  a  pre¬ 
diction  method  with  (s  =  10.495,  <  =  0.41,  K  =  .0168)(the  latter 
two  values  are  also  used  by  Cebeci  &  Smith,  1974);  however,  it 
was  determined  (after  the  prediction  method  was  compared  to  a 
number  of  data  sets,) that  the  velocity  profiles  were  not  well 
predicted,  particularly  in  flows  with  pressure  gradient.  The 
objective  of  this  study  is  to  investigate  a  systematic  way  of 
determining  any  trends  in  these  parameters  for  various  effects 
such  as  pressure  gradients  or  mainstream  turbulence, with  the 
ultimate  goal  of  obtaining  correlations  of  these  parameters 
for  use  in  a  prediction  method. 

This  will  be  carried  out  here  in  a  somewhat  different  way 
from  previous  investigations.  It  has  become  common  practice 
in  recent  times  to  determine  values  of  "universal  constants" 


in  turbulence  models  through  a  process  known  as  computer 
optimization;  this  procedure  is  often  associated  with  turbu¬ 
lence  models  that  have  a  sizeable  number  of  adjustable  constants 
(see,  for  example,  Murthy,  1977)  and  is  usually  carried  out 
as  follows.  One  data  set  or  a  number  of  data  sets  are  chosen, 
each  of  which  consists  of  a  number  of  measured  profiles  at 
various  downstream  locations  from  some  initial  point;  a  pre¬ 
diction  method  with  a  preassigned  set  of  constants  in  the 
turbulence  model  is  started  at  the  initial  data  station  and 
the  downstream  data  is  “predicted."  On  the  basis  of  compar¬ 
ison  with  either  the  integral  parameters  and/or  the  skin 
friction  coefficient,  a  decision  is  made  as  to  whether  adjust 
the  "constants"  in  the  turbulence  model  to  achieve  better 
"predictions";  the  prediction  method  is  then  used  to  obtain 
another  "prediction"  of  the  downstream  data.  This  iterative 
method  continues  until  some  optimal  set  of  constants  is 
obtained. 

One  undesireable  feature  of  this  scheme  is  that  the 
basis  for  adjusting  the  constants  is  usually  not  clearly 
defined  and  in  any  case  is  based  upon  comparisons  with  the 
gross  properties  such  as  the  integral  quantities  or  quan¬ 
tities  obtained  indirectly  from  experiment  such  as  the  skin 
friction.  The  details  of  how  the  iterative  process  is  carried 
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j  out  and  what  criteria  are  used  to  decide  on  the  optimal  set  of 

j  constants  are  normally  not  supplied  in  the  literature,  parti- 

'  cularly  when  the  number  of  data  sets  involved  are  large;  pre- 

i 

sumeably  the  choice  is  based  ultimately  on  some  type  of  inte- 
;  grated  subjective  average. 

i 

!  In  the  present  study,  the  trends  in  the  turbulence  model 

i 

1 

j  parameters  and  eventually  the  correlations  are  also  determined  by 

comparison  to  experimental  data  but  the  procedure  is  differ- 

I 

ent  and  precisely  defined.  First,  only  equilibrium  data  sets 
are  considered;  an  equilibrium  turbulent  boundary  layer,  by 
definition,  is  a  boundary  layer  in  which  the  mainstream  velo¬ 
city  varies  algebraically  with  streamwise  distance.  Such  a 
boundary  is  expected  (Clauser,  1956;  Mellor,  1972;  Fendell, 
1972;  Scharnhorst  et  al.,  1978)  to  successively  approach  a 
}  self-similar  flow  at  large  distances  from  wherever  the  boundary 

I  layer  flow  is  initiated.  For  this  reason,  equilibrium  boundary 

layers  have  historically  been  of  considerable  interest. 

It  is  worthwhile  to  remark,  however,  that,  in  order  for 

I 

self-similar  velocity  profiles  to  exist,  equilibrium  is  a 
necessary  but  by  no  means  a  sufficient  condition.  To  expand 
on  this  point,  consider  the  case  of  laminar  boundary  layers 

where  similar  solutions  satisfy  the  well  known  Falkner-Skan 

i 

I  equation.  The  existence  of  such  solutions  in  laminar  flows 
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has  been  discussed  by  Brown  &  Stewartson  (1965)  who  argue  that 
similar  solutions  may  be  expected  in  two  physical  situations. 
The  first  of  these  is  at  an  x  station  where  the  velocity 
profile  is  an  initiator  of  the  boundary-layer  flow  downstream; 
such  a  situation  occurs  physically  at  the  front  stagnation 
point  of  a  bluff  body  or  the  leading  edge  of  a  flat  plate,  for 
example.  In  these  cases,  the  Falkner-Skan  profile  gives  the 
proper  laminar  boundary- layer  solution  at  a  point  of  attach¬ 
ment  of  the  mainstream  and  provides  the  initial  condition  to 
initiate  a  boundary- layer  calculation  downstream.  For  tur¬ 
bulent  boundary  layer..,  there  appears  to  be  no  analogue  of 
this  physical  situation.  In  the  absence  of  mainstream  turbu¬ 
lence,  the  flow  at  a  point  of  attachment  of  the  mainstream  is 
observed  to  be  laminar  and  when  the  downstream  boundary- layer 
flow  is  turbulent,  there  is  a  transition  zone  in  between  the 
laminar  and  the  turbulent  flow.  Moreover,  a  wide  variety  of 
upstream  experimental  conditions  can  lead  to  transition  and  an 
eventual  fully-developed  turbulent  boundary  layer  downstream 
for  the  same  mainstream  velocity  distribution.  Consequently, 
there  would  appear  to  be  no  reason  to  expect  that  the  flow  in 
the  outer  region  of  an  equilibrium  boundary  layer  will  be  self¬ 
similar  at  the  initial  stations  of  a  fully  developed  turbulent 
flow. 
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For  laminar  boundary  layers,  the  other  case  discussed  by 
Brown  and  Stewartson  (1965)  is  that  where  a  similarity  solution 
becomes  what  may  be  described  as  a  terminator  of  a  more  general 
boundary- layer  flow  and  two  cases  of  this  behavior  are  con¬ 
sidered  by  Brown  and  Stewartson  (1965).  A  physical  case  where 
this  can  occur  is  at  a  point  of  detachment  of  the  inviscid  flow; 
examples  of  such  behavior  are  known  in  magneto- hydrodynamic  flow 
(Leibovich,  1967a, b)  and  in  rotating  flows  (Walker  &  Stewartson, 
1972)  at  the  rear  stagnation  point  of  symmetrical  and  two- 
dimensional  bluff  bodies.  Another  case  is  flat  plate  flow  and 
here  if  the  initial  velocity  profile  at  any  arbitrary  location 
on  the  plate  is  not  given  by  the  Blasius  solution,  then  the 
Blasius  profile  can  only  become  the  relevant  solution  at  an 
infinite  distance  downstream.  For  turbulent  boundary-layer 
flows,  an  analogous  type  of  situation  is  expected;  that  is, 
similarity  solutions  are  only  anticipated  as  terminators  and 
usually  at  large  distances  downstream  from  wherever  the  boun¬ 
dary  layer  is  initiated. 

In  practice,  one  would  expect  to  be  able  to  measure  tur¬ 
bulent  boundary-layer  profiles,  at  large  distances  downstream 
of  the  transition  zone,  which  become  arbitrarily  close  to 
being  self-similar;  however,  near  the  transition  zone,  there 
is  no  reason  to  expect  a  self-similar  behavior.  In  zero 
and  favorable  pressure  gradients,  measured  profiles  in  an 
equilibrium  flow  should  increasingly  approach  self-similarity 
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at  subsequent  data  stations  downstream.  In  adverse  pressure 
gradient  flows,  there  is  an  additional  difficulty  in  that 
reversed  flow  and  boundary  layer  separation  may  occur  before 
similarity  is  achieved. 

The  self-similar  flow  is  particularly  attractive  insofar 
as  turbulence  model  development  is  concerned  because  the  para- 

I 

bolic  partial  differential  equation  governing  the  turbulent 
boundary  layer  flow  becomes  an  ordinary  differential  equation 
of  the  boundary  value  type.  In  the  present  study,  the  ordinary 
differential  equation  associated  with  the  outer  layer  is  solved 
numerically  for  a  given  set  of  the  turbulence  model  parameters 
(k,K)  to  determine  an  outer  layer  velocity  profile;  this  outer 
layer  profile  is  matched  with  an  inner  layer  profile  containing 
the  inner  profile  parameters  (s,k)  and  a  composite  velocity 
profile  for  the  entire  boundary  layer  is  defined.  The  com¬ 
posite  profile  is  then  compared  directly  to  measured  experi¬ 
mental  velocity  profile  data  and  a  r  ot-mean-square  error  is 
defined  as  a  criterion  of  how  well  the  profile  represents  the 
data.  The  optimization  procedure  then  adjusts  one  or  more 
of  the  three  profile  parameters  (s,k,K)  until  a  "best  fit" 
to  each  profile  is  obtained.  Once  this  procedure  has  been 
carried  out  for  several  data  stations,  the  results  are  plotted 
to  determine  any  trends  in  the  profile  parameters;  the 
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objective  here  is  to  obtain  correlations  for  the  turbulence 
model  parameters  for  use  in  a  prediction  method.  In  the 
present  study,  this  procedure  is  carried  out  for  two  situations 
and  correlations  are  obtained  for  k  and  K  for  the  effects  of 
pressure  gradients  and  mainstream  turbulence.  Note  that  this 
type  of  procedure  is  only  applicable  to  measured  profiles  for 
which  the  mainstream  velocity  is  of  the  equilibrium  type  and 
where  the  data  is  at  locations  far  downstream  of  wherever  the 
turbulent  boundary  layer  was  initiated.  The  final  correla¬ 
tions  for  the  effects  of  pressure  gradients  and  mainstream  tur¬ 
bulence  lead  to  excellent  representations  of  the  velocity  pro¬ 
file  data.  In  principle,  the  models  developed  here  may  be 
used  in  a  prediction  method  for  non-equilibrium  flows. 

The  plan  of  this  report  is  as  follows.  In  §2,  the  basic 
equations  and  principle  results  of  the  asymptotic  theory  are 
sutimari zed .  In  §3  the  optimization  procedure  is  described 
and  in  §4  and  §5  the  results  of  the  method  for  pressure  gradients 
and  mainstream  turbulence  effects  respectively  are  given.  A 
description  and  test  cases  for  the  optimization  code  are  given 
in  Appendices  C  and  D.  Finally,  the  conclusions  of  the  study 
are  discussed  in  §6. 
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2.  TtiE  TURBULENT  BOUNDARY  LAYER  TIME  MEAN  EQUATIONS 


2.1  Basic  Equations 

The  turbulent  boundary  layer  equations  governing  two- 
dimensional,  incompressible  time-mean  flow  are  the  mass  continu¬ 
ity  equation  and  the  Reynolds  equation  for  streamwise  momentum: 


M  +  0 

3x  ay  ’ 


n  lu  - 


-  au  . 
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(2.1) 

(2.2) 


In  these  equations,  (x,y)  represent  Cartesian  coordinates  with 

corresponding  mean  velocities  (u,v).  All  lengths  and  velocity 

components  have  been  made  dimensionless  with  respect  to  L  and  Uq, 

a  reference  length  and  velocity  respectively.  The  Reynolds  num- 
UqL 

ber  is  defined  as  and  is  assimed  to  be  large;  here  v  is  the 
kinematic  viscosity.  The  momentum  equation  (2.2)  contains  an 
additive  stress  term  o  for  turbulent  flow  which  is  related  to 
the  usual  turbulent  shear  stress  by 


a  = 


-uV 


(2.3) 


In  equation  (2.3)  the  primed  quantities  are  instantaneous  velo¬ 
city  fluctuations  about  the  corresponding  mean  values  and  the 
over  bar  signifies  a  long  time  average.  The  boundary  conditions 


associated  with  the  equations  (2.1)  and  (2.2)  are 


u(x,0)  =  v(x,0)  =  0  , 
u(x,y)  ^  U„(x)  as  y  . 


(2.4) 


In  addition  to  the  above  boundary  conditions,  the  turbulent 
shear  stress  term  must  be  chosen  to  satisfy 


a(x,0)  =  0,  a(x,y)  -»■  0  as  y  -►  »  .  (2.5) 

To  completely  specify  the  problem,  an  initial  velocity  profile 

u(0,y)  =  f(y)  for  0  £  y  <  »  ,  (2,6) 

must  be  given  at  some  x-station  in  the  fully  turbulent  region 
of  the  flow  which  is  denoted  here  by  x»0.  Note  that  the  problem  is 
indeterminate  until  a  model  for  the  Reynolds  shear  stress  term 
is  specified;  however,  it  is  possible  through  asymptotic  analy¬ 
sis  of  this  problem  to  provide  useful  information  about  the 
velocity  components  and  shear  stress  without  the  introduction 
of  any  specific  functional  form  for  the  Reynolds  stress  term. 
Before  this  is  carried  out,  it  is  desirable  to  define  an  impor¬ 
tant  physical  parameter,  u*,  that  will  play  an  important  role 
in  the  asymptotic  analysis.  The  dimensional  wall  shear  stress 
is  given  by 
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(2.7) 


and  is  used  to  define  the  dimensional  skin  friction  velocity 


“t  *  *  (2.8) 

where  and  represent  the  dynamic  viscosity  and  density  at 
the  wall  respectively.  A  dimensionless  skin  friction  velocity 
can  be  defined  as 


* 
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u 

uX(x) 


(2.9) 


Note  that  equation  (2.7)  can  be  rewritten  in  an  equivalent  form 
as 


1  au 


(2.10) 


The  asymptotic  structure  of  equation  (2.2)  in  the  limit  of  large 
Reynolds  number  will  now  be  considered. 


2.2  Asymptotic  Structure  of  the  Time-Mean  Equations 

It  is  well  known  that  the  turbulent  boundary  layer  is  a 
composite  double  layer  consisting  of  a  thin  inner  layer  adjacent 
to  the  wall  and  a  relatively  thick  outer  layer.  A  number  of 
authors,  including  Fendell  (1972)  and  Mellor  (1972),  have 
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considered  the  asymptotic  structure  of  the  time-mean  momentum 
equations  for  a  constant  property  incompressible  flow;  in  these 
studies,  asymptotic  methods  are  used  to  isolate  the  leading  terms 
in  an  expansion  for  large  Reynolds  number  for  the  flow  quantities 
in  both  the  inner  and  outer  region  of  the  turbulent  boundary 
layer.  These  expansions  are  then  matched  in  a  manner  which  is 
consistent  mathematically  and  which  is  also  consistent  with 
experimental  measurements  of  the  time-mean  flow  quantities  for 
both  regions.  Fendell  (1972)  also  examines  the  conditions 
necessary  for  self-similar  flow  in  the  outer  region.  The  asymp¬ 
totic  results  that  follow  summarize  the  work  of  Fendell  (1972) 
and  are  also  discussed  by  Scharnhorst  et  al  (1978). 

In  the  outer  region  of  the  turbulent  boundary  layer,  the 
velocity  profile  may  be  written  as  a  small  perturbation  about 
the  mainstream  value  as  Re  -*■  *»,  according  to, 

3F, 

u  =  U«(x){l+u*(x;Re)  (x,n)  +  (2.11) 

Here  u*  is  the  dimensionless  skin  friction  velocity  defined  by 
equation  (2.9)  and  u*-*0  as  Re^.  In  the  outer  region,  the  tur¬ 
bulence  term  o  may  be  written  to  leading  order  as 

a  ■  U„2(x){u*2(x;Re)Ei(n.x)+  ...}  .  (2.12) 

Here  n  is  the  scaled  normal  coordinate  for  the  outer  region 


defined  by, 


n  *  y/Ag 


(2.13) 


The  dimensionless  outer  region  length  scale  is  proportional 
to  the  boundary  layer  thickness  and  is  0(u*);  a  convenient 
choice  for  a^  is  made  in  section  2.4;  It  is  of  interest  to  note 
that  the  shear  stress  term  given  in  (2.12)  is  equivalent  to  the 
result  that  u'v'  is  0(u2^)  which,  in  general,  is  confirmed  by 
experimental  measurement.  Vihen  the  velocity  profile  expansion 
(2.11)  and  the  turbulence  term  expansion  (2.12)  is  substituted 
into  (2.2)  and  terms  quadratic  in  the  perturbation  u*  are  neg¬ 
lected,  an  equation  for  the  velocity  defect  function  u-j  =  3F^/3n 
is  obtained  according  to 


311  _  A^ 

IT'  7^ 


(u*U»2)- 


(U«Ap)'n 


U*"  3ti3X  ■ 


(2.14) 


Here  a  prime  denotes  differentiation  with  respect  to  x. 

This  equation  is  subject  to  the  boundary  conditions 

Fl(x,n)  ->•  0  as  n  -*•  0  ,  (2.15) 

and 

3F, 

(x,n)  ^  0  as  n  *  .  (2.16) 
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In  the  inner  sometimes  referred  to  as  the  wall  layer 

the  velocity  profile  and  turbulence  stress  term  expansions  may 
be  written  to  leading  order  according  to 

u  =  u  U«,(x)  ^  (x.y'*’)  +  ...  ,  (2.17) 

and 

a  =  U«2(x)u*(j^(x,y+)  +  ...  ,  (2.18) 

where 

.  u  L 

y‘*‘  =  ^y  =  Re  u*U«>(x)y  .  (2.19) 

is  the  inner  region  variable.  The  matching  of  the  leading  order 
terms  in  the  asymptotic  expansions  (2.11)  and  (2.17)  and  (2.12) 
and  (2.18)  occur  in  the  limits  ty+0  and  y^-^*  for  the  inner  and 
outer  layers  respectively.  A  self  consistent  mathematical 
structure  which  is  compatible  with  experimental  measurement  may 
be  obtained  if  both  velocity  profiles  merge  smoothly  with  a 
logarithmic  profile  behavior  according  to 

u-,  -  log  n  +  Cq(x),  as  n  0  ,  (2.20) 

and 

u"*"  -  log  y"*”  +  (x),  as  y"*"-^.  (2.21) 

Here  k,  C  and  are  in  general  functions  of  x.  Although  < 


■Tl 


is  analogous  to  the  von  Karman  “constant"  which  is  normally 
assumed  to  have  a  value  of  about  0.41,  in  general  k  could 
depend  on  local  flow  conditions.  It  is  worthwhile  to  note  that 
the  conditions  for  the  turbulence  terms  also  follow  from  the 
analysis  (Fendell,  1972)  and  are 

1  as  ri  -►  0  ,  (2.22) 

and 

o-j  ->•  1  as  y"*"  -►  «  .  (2.23) 

First  order  matching  of  the  inner  and  outer  asymptotic 
expansions  leads  to  the  velocity  match  condition  given  by. 


W  “  ^  Re  u*U»(x)}  +  .  (2.24) 


This  skin  friction  relation  connects  u*  and  the  outer  region 
length  scale  The  match  condition  to  leading  order  as  Re  » 
implies  that 


u 


* 


<(xl 

log  Re 


(2.25) 


and  since  is  0(u*)  the  turbulent  boundary  layer  thickness 
approaches  zero  as  the  inverse  of  a  logarithm  in  the  high 
Reynolds  number  limit. 

A  composite  velocity  profile  valid  to  leading  order  across 
the  entire  boundary  layer  can  be  formulated.  Van  Dyke  (1975) 
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defines  a  composite  expansion  according  to 


[composite  expansion]  =  [inner  expansion] 

+  [outer  expansion]  -  [common  terms].  (2.26) 

The  common  terms  represent  the  behavior  of  the  inner  and  outer 
expansions  in  the  matching  region.  A  composite  streamwise 
velocity  profile  is  defined  here  according  to, 

u  =  u*U„[u+  +  ^  log  n  +  }J.  (2.27) 

Alternatively,  an  equivalent  expression  could  be  composed  by 
using  equation  (2.24)  for  the  common  terms  in  terms  of  the  inner 
region  variable  y+. 

The  question  of  determining  good  model  profile  approxima¬ 
tions  for  the  inner  profile  U+  and  the  outer  profile  will 
now  be  considered. 

2.3  Similarity  in  the  Outer  Region 

The  model  profile  that  will  be  used  for  the  outer  region 
flow  is  a  self-similar  profile  and  for  this  reason  it  is  of 
interest  to  examine  the  conditions  necessary  for  self-similar 
flow  in  the  outer  region.  The  terms  self-similar  and  equilibrium 
are  often  used  interchangeably  in  the  literature  but  it  is  important  to 
make  a  distinction  here.  The  term  equil  ibrium  is  understood  to  apply 
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to  a  turbulent  boundary  layer  for  which  the  mainstream  velocity 
behaves  as  Uoo(x)  ~  x“  or  Uco(x)  ~  e^’^  where  a  and  g  are  constants; 
on  the  other  hand,  self  similarity  of  the  outer  layer  velocity 
profile  can  only  occur  when  the  boundary  layer  is  exposed  to 
the  equilibrium  outer  velocity  contributions  for  large  stream- 
wise  distances.  Consequently,  equilibrium  is  a  necessary  but  not 
sufficient  condition  for  the  existence  of  self-similar  profiles. 

The  necessary  conditions  for  similarity  follow  from  equa¬ 
tion  (2.17)  for  and  are  that. 


=  -2g  =  constant  (2.28) 

and 

-;Jt-  (a*  U»)'  =  a  =  constant.  (2.29) 


To  assess  the  magnitude  of  the  ratio  u*7u*  which  appears  in 
equation  (2.28),  the  match  condition,  given  by  equation  (2.24), 
is  differentiated  with  respect  to  x  to  obtain, 

A  *  U  * 

■-logU^Re  u*U.(x))  +  7  (-f*  +  '-C„' • 

(2.30) 


-K 

7^ 


This  equation  can  be  simplified  by  using  the  match  condition  (2.24) 
to  eliminate  the  first  term  on  the  right  hand  side  and  this  pro¬ 
cedure  yields 
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u*' 

u* 


+  k'(C-Co)+<(C.-C^)'}. 

(2.31) 


Since  k  must  approach  a  constant  for  a  self-similar  flow,  the 
order  of  magnitude  of  the  ratio  u*'/u*  follows  from  equation 
(2.31)  and 

^“0(u*).  (2.32) 


Terms  containing  this  ratio  may  be  neglected  to  leading  order 
in  (2.28)  and  (2.29)  which  become 

U  ' 

q-if-=-6  (2.33) 

^00 

and 


q‘ 


q  =  a 


(2.34) 


respectively  where  q  *  a^/u*.  These  two  equations  are  combined 
to  give 

q'  =  a  +  8  .  (2.35) 


As  a  result,  there  are  only  two  types  of  mainstream  velocity 
distributions  which  can  lead  to  self-similar  solutions  in  the 
outer  layer  and  these  are: 

u«(x)  =  D^(x-x^j)“^  for  a  +  8  0  (2.36) 
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where 


-1  =  •  <1  '  #  '  (■»+6)(x-x<,)  (2.37) 


.  ^0  _ 


where 


Ua.(x)  =  D2e  for  a  +  B  =  0  (2.38) 


=  ^  .  q  =  717  =  •  (2.39) 


Here  D.|,  D2,  and  are  all  arbitrary  constants. 

Equilibrium  flows  have  been  examined  experimentally  for  a 
number  of  years  and  an  equilibrium  boundary  layer  has  been 
defined  experimentally  as  a  flow  in  which  the  dimensionless  velo¬ 
city  defect  (Ug-  u*)/u^  expressed  as  a  function  of  y/S  becomes 
close  to  being  invariant  with  downstream  distance.  After  careful 
experimentation,  Clauser  (1954,1956)  concluded  that  the  criterion 
for  equilibrium  was  a  constant  value  of  the  parameter  e  which 


is  defined  as 


A* 

^c  iF~  *^e  "djT 

T 


(2.40) 


Here  6  represents  the  dimensionless  displacement  thickness 


defines  as 


so 


(2.41) 


From  the  definition  of  6*  and  equations  (2.11)  and  (2.13)  it 
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follows  that. 


5*  =  -AqU*F^(x,»)  +  . . . 


(2.42) 


As  self-similarity  conditions  are  approached,  the  x  dependence 
must  vanish  and  Fi(x,«)  must  approach  a  constant  value,  say  F-iC"). 
Thus  (2.40)  becomes 


& 


c 


Aq  F-|(“)  du„ 

u*Uoo  djT 


=  -3  F^(oo)  . 


(2.43) 


Note  that  a  constant  value  of  3  implies  a  constant  value  of  3_ 
and  therefore  Clauser's  (1954,1956)  experimental  results  are 
consistent  with  theory. 


2.4  The  Eddy  Viscosity  Model  for  the  Outer  Region 

In  order  to  obtain  a  solution  of  the  outer  layer  equations, 
a  constitutive  relation  for  u'v*  in  the  outer  layer  must  be 
specified  and  for  the  reasons  discussed  in  section  1,  a  simple 
eddy  viscosity  model  will  be  used  here.  It  is  worthwhile  to 
note  that  although  it  is  customarily  assumed  that  some  functional 
relationship  exists  between  u 'v'  and  the  mean  velocity  profile, 
no  such  relationship  has  been  demonstrated  either  experimentally 
or  theoretically;  consequently,  the  eddy  viscosity  hypothesis 
should  be  regarded  at  present  as  a  convenient  approximation 
which  is  expected  to  be  supplanted  in  the  future  by  constitutive 
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models  based  on  the  observed  coherent  structure  of  the  turbulent 
boundary  layer. 

For  the  outer  region,  the  Reynolds  stress  is  defined  as 


(2.44) 


and  from  equations  (2.12)  and  (2.13),  the  stress  function  z 


becomes 


3U, 


'1  ~  AqU^Uo„(x)Lu* 


1 


(2.45) 


where  depends  upon  both  x  and  n.  The  functional  form  for 
is  the  same  type  as  that  used  in  the  Cebeci-Smith  (1969)  and 
Mellor  and  Herring  (1969)  prediction  methods.  In  this  model, 
Eq  is  selected  to  approach  a  constant  value  for  fixed  x  near 
the  outer  edge  of  the  outer  layer.  Thus, 


eo  K  Ug(x)6*L  (2.46) 

where  6*  is  the  dimensionless  displacement  thickness  and  K  is 
an  empirically  determined  constant.  The  value  of  K  differs 
slightly  according  to  the  model;  Cebeci  and  Smith  (1969)  take 
K  =  0.0168  while  Mellor  and  Herring  (1968)  use  K  =  0.016.  In 
the  present  study,  a  universal  value  for  K  is  not  assumed;  rather 
one  objective  is  to  determine  if  this  parameter  depends  in  any 
way  on  the  pressure  gradient. 
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Near  the  inner  edge  of  the  outer  layer,  the  eddy  viscosity 


must  behave  linearly  to  satisfy  the  matching  condition  for  E-j 
given  in  equation  (2,22);  it  is  convenient  to  define  a  dimension 
less  eddy  viscosity  function  E(n)  according  to 


"  6%u!(x)L  * 
whereupon  equation  (2.45)  becomes 


(2.47) 


3U, 


'1  -  iT^  — 


(2.48) 

The  function  £(n)  must  assiine  the  following  limiting  values: 

(2.49) 


u*A, 


and 


e  <  n  (  Tjr  )  as  n  -»■  0  , 

0 


e  ->■  K  as  n  ->■ 


(2.50) 


A  particularly  convenient  choice  for  the  outer  region 
length  scale  is 


A 


0 


(2.51) 


and  it  follows  that  F^(x,«»)=-1  from  equation  (2.42).  Detailed  reasons 
behind  this  choice  for  Aq  are  discussed  by  Scharnhorst  et  al. 
(1978)  and  Scharnhorst  (1978)  and  Weigand  (1978).  Equations 
(2.45)  and  (2.49)  may  now  be  rewritten  as 


(2.52) 


and 


3U^ 

=1  ■  1^" 


e  -»•  icn  as  n  . 


(2.53) 


In  general,  the  functional  form  of  e(n)  must  be  specified.  The 
model  of  Mel  lor  and  Gibson  (1966)  and  Mel  lor  and  Herring  (1969) 
use  a  simple  form  of  two  straight  lines  defined  by 

e(n)  =  e_(n)  =  <f  for  n  mi  .  (2.54) 

[  icn  for  n  <  Hj  , 

where  n-j  =  K/k.  An  awkward  feature  of  this  model  lies  in  the  dis¬ 

continuous  derivative  at  n  =  n-i  which  may  be  expected  to  give 
rise  to  difficulties  in  a  numerical  solution  of  the  outer  layer 
equation;  this  problem  is  handled  by  Cebeci  and  Smith  (1974) 
by  using  an  artificial  smoothing  of  the  model  in  equation  (2.54). 
An  alternative  functional  form  for  e(ri)  must  reflect  the  linear 
behavior  near  n=0  in  equation  (2.53)  and  the  limiting  constant 
value  in  equation  (2.50)  for  large  n;  moreover,  such  a  function 
should  be  monotonically  increasing  with  n.  A  rather  complex 
function  meeting  these  requirements  was  assumed  by  Scharnhorst 
et  al.  (1977)  who,  in  addition,  determined  that  it  appeared  to 
be  desireable  for  e(n)  to  approach  the  linear  behavior  icn  expo¬ 
nentially  quickly  for  small  values  of  n.  It  also  appeared 
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important  for  e(n)  to  approach  the  upper  bound  K  relatively 
quickly  for  large  n. 

A  particularly  simple  form  of  a  function  that  satisfies 
all  the  desired  features  of  e(n)  is 
,  ,  -C/n*^  1/N 

e(n)  =  Kn{l-e  }  (2.55) 

which  is  plotted  in  Figure  2.1  for  integral  values  of  N  =  1,  2, 
and  4.  Here,  the  term  C  is  given  by 

C  =  (f)  .  (2.56) 

This  monotonically  increasing  function  (2.55)  meets  the  require¬ 
ments  as  n  -*■  0  , 

e(n)  -*■  Kn  .  (2.57) 

and  as  n  -►  “  » 

E(n)  K{1  +  0(4-)}  .  (2.58) 

Note  that  this  function  for  N=4  is  virtually  identical 
to  the  form  used  by  Schamhorst  et  al.  (1977)  but  is  a  much 
simpler  form. 

The  outer  region  eddy  viscosity  function  given  by  equation 
(2.55)  was  used  throughout  this  study  with  N=4. 
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Viscosity 


With  the  outer  layer  stress  function  defined  by  equation 
(2.52),  it  may  be  shown  that  the  outer  layer  equation  (2.14) 


becomes 

3n  3n  ^  u*U»  ^  3ti  u*^U«^  ^1  u*’  3x  ’ 

(2.59) 

where  *  3F^/3n.  For  a  self-similar  flow,  equation  (2.59) 
reduces  to  the  ordinary  differential  equation 


H  dU, 

^  “d^f^  ^ 


.6*' 


dU, 


‘dRT  ^ 


U,  =  0 


(2.60) 


where  6^  is  the  Clauser  pressure  gradient  parameter  obtained 
from  equation  (2.40)  as  =  -6*U»7(u*^Uco).  Integration  of 
equation  (2.60)  across  the  boundary  layer  yields  the  relation 

^  =  1  +  3b^  ,  (2.61) 

which  can  be  substituted  back  into  (2.6i,  to  obtain 
H  ‘lU,  dU, 

^  (c(n)  -3^)  +  (H2s^)n  -jJ-  +  26^  U,  •  0.  (2.62) 

To  find  the  velocity  profile  in  the  outer  region,  equation 
(2.62)  must  be  solved  subject  to  the  boundary  conditions , 
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3s  n  0  » 


(2,63) 


U,  -  1  log  .  ♦  C„ 

U.|  ->■  0  as  n  .  (2.64) 

Note  that  the  constant  is  unknown.  This  problem  was  solved 
by  a  combination  of  a  series  solution  near  ii=0  and  a  fully  numer¬ 
ical  solution  for  n>0;  the  procedure  is  described  in  Appendix  A. 

Note  that  there  are  two  parameters  associated  with  equation 
(2.62)  and  the  boundary  condition  (2.63);  these  are  the  eddy 
viscosity  parameters  K  and  the  von  Karman  constant  k  which  appears 
in  the  eddy  viscosity  and  in  the  boundary  condition.  In  addi¬ 
tion  to  the  physical  boundary  layer  quantities,  the  skin  friction 
u^  and  the  displacement  thickness  are  contained  in  the  parameter 
Bj.  which  is  defined  in  equation  (2.40)  and  which  appears  in  equa¬ 
tion  (2.62), 


2.6  Wall  Layer  Model 

The  wall  layer  model  used  in  this  study  incorporates  what 
are  believed  to  be  the  important  features  of  the  time  dependent 
flow  in  the  wall  layer.  A  complete  discussion  of  the  ideas  that 
develop  the  unsteady  wall  layer  model  is  given  in  Walker,  Abbott 
and  Scharnhorst  (1976),  Walker  and  Abbott  (1977),  and  Walker 
and  Scharnhorst  (1977).  This  "unsteacty  wall  layer  model"  has 
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been  extensively  compared  to  experimental  data  by  Scharnhorst 
et  al .  (1977)  who  demonstrate  that  the  model  closely  represents 
measured  velocity  profile  data  in  the  wall  layer  even  in  flows 
with  pressure  gradient;  moreover,  the  representation  is  con¬ 
siderably  better  than  that  obtained  by  a  conventional  Van  Driest 
(1956)  type  of  model  and,  finally,  the  model  is  simpler  to  use. 
For  these  reasons,  this  "unsteady  wall  layer  model"  was  used 
exclusively  in  the  present  study. 

The  profile  given  by  Walker  and  Scharnhorst  (1977)  is 


V 


=  [1  +  ^][R(S2,to+)Q(H)  +  Z(H)  +  P(S2,to+)W(H)] 
t-+ 

-  -jr  CR(0,t/)QCHQ)  +  Z(Hj,)  +  P(0,to+)W(HQ)],  (2.65) 


where 


H  * 


JL 


2/52+t^ 


(2.66) 


and 


R(S2,to+)  =  -  log  2}  +  ]-  p+(S2+2to+) 


+  ^  log  (S2+to+)  ,  (2.67) 

Q(H)  =  (2H2+1)  erf  H  +  -?-  H  e'”^  ,  (2.68) 

/x 

Z(H)  =  [(2H2+1)5(H)  +  H5'(H)-  (6H2+l)erf  H  -  |h  e"”^], 

C2.69) 


(2.70) 


P(S^t/)  =  -  I  P+(S2+t^+)  . 

W(H)  *  [HH3H2+  |]erf  H  -  —  [H2+  |]e“^^-3H2  .  (2.71) 

/iT 

In  the  equations  above,  the  H  function  is  defined  as 

h  X  y 

5(n)  =  I  \  \ 

0  0  0 

A  detailed  description  of  the  h  function  and  its  properties  is 
given  in  Appendix  B.  For  large  y+  in  the  inner  region,  the 
asymptotic  form  of  the  profile  (2.65)  is 

U*"  -  1  log  y+  +  C.  .  (2.73) 

Note  that  the  model  contains  k  and  as  parameters  in  addition 
to  tg''’  and  S;  these  last  two  parameters  are  associated  with 
the  physics  of  the  wall  layer  time  dependent  flow  and  are  dis¬ 
cussed  in  detail  by  Walker  and  Scharnhorst  (1977)  and  Scharnhorst 
(1978).  The  mean  velocity  profile  of  (2.65)  must  satisfy 

(2.74) 

(2.75) 
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3U+ 

W 


y+=0 


=  1 


and  the  wall  compatibility  condition 


93U+ 

3y+^ 


ly+sQ 


=  0 


^  ] 

I 

i 

which  leads  to  the  auxiliary  relations 

(S2+t„+)i[R(S2,t„*)-  1  +  PCS2,t/)] 

-  ^  +  P(0,t„-*)]  =  $2  .  (2.76) 

and 

(S2+to‘')‘^R(S2.t/)+  3P(S2.tQ+)] 

,  ~i  .  . 

-(t^  )  [RCO.t^  )  +  3P(0.tjj+)]  =  0  ,  (2.77) 

respectively.  Note  that  equations  (2.76)  and  (2.77)  are  two 
relations  for  the  four  parameters  ic,  C^-,  S,  and  t^^  and  conse¬ 
quently  only  two  of  these  parameters  are  independent.  In  the 
data  comparisons  that  are  carried  out  here,  k  and  S  were  gener¬ 
ally  assumed  to  be  variables  that  were  adjusted  to  obtain  the 
best  fit  to  the  data;  thus,  at  any  stage  in  the  optimization 
procedure  for  specified  values  of  k  and  S,  equations  (2.76) 
and  (2.77)  were  solved  for  to'*’  and  C.j. 

2.7  Summary 

A  composite  similarity  velocity  profile  valid  to  leading 
order  for  the  entire  boundary  layer  has  been  developed  in  this 
section  according  to 

I 

i 

u  *  u*U«p  +  U,  -  log  n  +  y]  .  (2,78)  j 
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Here  U+  is  the  unsteady  wall  layer  model  given  in  section  2,6; 
this  model  is  an  analytic  expression  given  by  equation  (2.65) 
which  contains  the  independent  parameters  k  and  S.  The  outer 
region  profile  U^(n)  must  be  obtained  as  a  numerical  solution  of 
differential  equation  (2.62)  developed  in  section  2.5;  this 
numerical  solution  will  implicitly  contain  the  parameters  k 
and  K  which  appear  in  the  eddy  viscosity  model.  Thus,  the  com¬ 
posite  profile  contains  the  three  independent  model  parameters 
S,  K  and  K  that  may  be  adjusted  in  a  computer  optimization  rou¬ 
tine  to  obtain  a  best  fit  with  experimental  velocity  profile  data. 
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3.  THE  OPTIMIZATION  TECHNIQUE 


3.1  Background 

The  present  turbulence  velocity  profile  model  contains 
three  adjustable  model  parameters  S,  k,  and  K  that  affect  the 
basic  shape  of  the  analytical  profile.  This  section  addresses 
the  development  of  a  systematic  method  to  determine  the  optimum 
values  of  the  model  parameters  that  minimize  the  difference 
between  the  analytical  model  profile  and  experimental  data. 

The  method  described  here  can,  in  principle,  be  applied  to  any 
test  data  to  attempt  to  determine  any  trends  in  the  data  with 
different  physical  effects;  in  particular,  the  procedure  is 
applied  here  to  data  for  several  pressure  gradient  flows  and 
mainstream  turbulence  levels. 

Before  examining  the  optimization  method,  a  basis  for 
determining  a  best  fit  must  be  defined.  For  this  study  a  root- 
mean-square  error  e  was  selected  according  to 


r2  = 


1  N  f 
^  n=l  [ 


'OATA'-'n 


iyj  -  u 


ANALYTICAL 


U. 


where 

N  =  number  of  data  points 

“data  "  experimentally  measured  velocity  at  y^^ 

“analytical  ^  analytical  velocity  profile  at  y^ 
Ug  =  freestream  velocity  . 


(3.1) 


The  best  analytical  fit  to  a  given  experimental  profile  is 
defined  as  that  set  of  model  parameter  values  which  minimizes 
e  as  defined  by  (3.1).  Previous  prediction  techniques  such 
as  those  presented  in  Scharnhorst  et  al.  (1977)  and  those  used 
in  the  Stanford  Conference  (Coles  &  Hirst,  1969)  utilize  the 
standard  error  criterion  of  (3.1)  as  an  objective  basis  of 
comparison.  The  particular  optimization  method  used  in  this 
study  is  presented  in  the  next  subsection. 

3.2  The  Optimization  Method  -  Direct  Search 

A  direct  search  minimization  procedure  was  used  to  optimize 
the  composite  velocity  profile.  This  technique  is  simple  to  use 
since  it  only  requires  values  of  the  objective  function  and  not 
gradients  to  carry  out  the  optimization  search.  The  nature  of  the 
problem  is  such  that  gradients  of  the  objective  function  cannot  be 
computed  analytically  and  can  only  readily  be  evaluated  by 
numerical  differentiation.  Although  the  direct  search  pro¬ 
cedure  becomes  very  time  consuming  when  the  number  of  optimiza¬ 
tion  variables  is  large,  it  was  used  here  with  good  success 
since  the  maximum  nunber  of  variables  is  three.  The  basic 
procedure  in  the  direct  search  minimization  was  to  vary  one 
model  parameter  such  that  the  least  squares  error  objective 
function  would  continuously  decrease  while  holding  the  other 
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parameters  constant.  Each  model  parameter  is  varied  in  turn 
so  that  at  any  stage  in  the  optimization  procedure,  a  one¬ 
dimensional  search  is  made  for  the  minimum  in  that  direction. 

The  method  used  here  is  somewhat  similar  to  the  direct  search 
method  of  Hooke  and  Jeeves  which  is  described  in  Hinmeblau  (1972) 
p.  142.  The  specific  logic  involved  is  as  follows. 

The  optimization  schane  starts  with  initial  values  for 
the  model  parameters  which  must  be  provided  to  the  subroutine 
as  well  as  a  vector  of  initial  incremental  step  sizes  for  the 
parameters.  This  initial  location  is  established  as  a  base 
point.  To  initiate  a  search,  the  objective  function  f°(x)  is 
evaluated  at  this  base  point  and  one  parameter  is  then  incre¬ 
mented  by  the  specified  step  size  Ax-j.  Suppose  first  that  the 
objective  function  decreases.  The  parameter  is  incremented 
continually  according  to  +  Ax^  and  the  objec¬ 
tive  function  is  computed.  This  process 

continues  for  i  =  0,  1,  2,  3,  ...  as  long  as  f^^^  continually 
decreases.  Eventually,  at  a  certain  step,  say  step  k,  the 
objective  function  will  increase  and  At  this 

stage,  a  local  minimum  in  the  x-j  direction  has  been  bracketed 
in  the  range  further  refine 

the  location  of  this  minimum,  a  quadratic  interpolation  poly¬ 
nomial,  f  =  Ax-j^  +  Bx-j  +  C  is  used,  where  A,  B,  and  C  are 
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determined  by  solving  the  set  of  equations: 


where 


+  C 


i  =  k-1,  k,  and  k+1 


(3.2) 


This  interpolation  polynomial  produces  a  local  minimum  in  the 
X"!  direction  at  x^*  =  -B/2A.  In  certain  cases, this  inter¬ 
polation  scheme  may  produce  a  value  of  x-j*  which  is  not  close 

fkl 

to  the  midpoint  of  the  interval  x-j'  ‘  or  an  objective  function 

/*)  /bi 

which  f'  '  is  not  close  to  f'  When  this  situation  occurs, 

then  three  of  the  values  x^*,  and 

which  have  the  lowest  objective  function  are  relabeled  as 

and  and  another  quadratic  interpolation 

polynomial  is  fit  through  these  values  to  obtain  a  new  minimum 
x-i*.  The  interpolation  scheme  is  repeated  until  a  specified 
convergence  criteria  is  met  between  the  values  of  x-j*  and 

or  f^^^  and  After  convergence  in  the  x-j  parameter, 

the  new  value  of  x-j  is  retained  as  the  base  point;  the  next 
parameter  x^  is  incremented  by  the  specified  step  size  AX2, 
and  the  search  procedure  repeats  until  all  the  independent 
parameters  have  been  changed. 

In  the  second  place  suppose  that  an  increase  in  the  objec¬ 
tive  function  is  realized  for  the  initial  step  of  a  model 
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parameter  (for  example,  x)  from  its  base  point.  When  this  situa¬ 
tion  occurs,  the  search  direction  is  inmediately  changed  to  the 
negative  direction  and  the  program  procedure  continues  as  pre¬ 
viously  stated.  In  the  event  that  the  objective  function  increases 
for  the  first  step  in  the  negative  direction  as  well,  the  program 
enters  into  the  quadratic  interpolating  routine  since  the  local 
minimum  has  been  bracketed  in  the  range  Xi-Ax.|  £  x-j  £  x-j  +  ax-j. 
After  the  local  minimum  has  been  refined,  a  one-dimensional 
search  is  carried  out  for  the  next  variable. 

After  a  one-dimensional  search  for  each  variable  has  been 
carried  out  by  the  above  procedure,  a  new  base  point  is  estab¬ 
lished  and  a  convergence  test  is  performed  in  which  the  optimiza¬ 
tion  function  values  for  the  last  two  base  points  are  compared. 

If  these  two  values  differ  by  an  amount  less  than  a  specified 
criterion,  a  minimum  has  been  found  and  the  program  terminates. 

If  the  convergence  criterion  is  not  met,  the  step  sizes  AXi 
of  the  search  are  reduced  ten  percent  and  the  search  procedure 
starts  anew  from  the  current  base  point. 

The  direct  search  program  contains  two  types  of  error 
flags.  The  first  type  of  error  flag  is  encountered  when  the 
number  of  combined  step  and  interpolation  iterations  reaches 
a  specified  maximum  value.  This  error  flag  then  terminates 
the  optimization  routine  and  returns  to  the  prediction  code 


with  the  latest  values  for  the  optimized  parameters.  The 
second  type  of  error  flag  occurs  in  the  quadratic  interpolation 
scheme  when  the  value  of  the  A  coefficient  of  the  quadratic 
polynomial  is  identically  zero.  This  flag  prevents  a  compu¬ 
tation  error  in  the  interpolation  process  and  terminates  the 
direct  search  program.  One  variable  returned  by  the  subroutine 
is  the  error  variable  lER;  values  lER*  1  and  lER  =  2  indicate 
errors  of  the  first  and  second  type  respectively.  A  value 
lER  =  0  indicates  a  successful  search  has  been  completed. 
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4.  TURBULENT  BOUNDARY  LAYER  PROFILE  REPRESENTATION 
4.1  Introduction 

In  order  to  obtain  a  velocity  profile  representation  for  a 
given  set  of  velocity  profile  data,  one  or  more  of  the  composite 
profile  parameters  S,  k,  and  K  can  be  alterect  in  several  combina¬ 
tions.  Each  such  combination  defines  a  different  method  of  pro¬ 
file  representation;  if  a  parameter  is  not  optimized,  it  must  be 
assigned  a  universal  value.  In  the  study  of  Scharnhorst  (1978), 
one  method  considered  was  to  take  ic=0.41  (the  generally  accepted 
universal  value  for  the  von  Karman  constant)  and  $=10.4965;  this 
value  of  S  is  the  value  which  with  ic=0.41  produces  a  value  of 
C^-  =  5.1  in  equations  (2.76)  and  (2.77);  again  C^-  =  5.1  is  a 
generally  accepted  value,  Scharnhorst  (1978)  then  carried  out  a 
one  parameter  optimization  on  K  over  a  wide  range  of  pressure 
gradient  data;  the  results  of  this  optimization  were  not  encourag¬ 
ing.  The  principal  difficulty  was  that  the  analytical  profile 
tended  to  skew  through  the  data  in  the  logarithmic  zone.  This 
difficulty  was  noted  by  Scharnhorst  (1978)  and  is  illustrated  in 
Figure  4.1;  in  this  figure  the  results  of  a  one  parameter  opti¬ 
mization  on  K  (with  <=0.41,  S»10.4965)  are  illustrated  for  three 
stations  of  the  data  of  Anersen  et  al  (1972).  These  stations 
are  labeled  8109,  8209,  and  8309  and  are  the  last  measured  data 
stations  in  equilibrium  flows  for  zero,  moderate  and  strong 
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FIGURE  4.1  Velocity  profile  comparisons  for  one  parameter 
optimization  on  K  (with  k»0.41,  8=10.4965)  for 
three  data  stations  of  Andersen  et  al.  (1972). 
Profiles  are  for  a  zero,  mild  adverse  and  strong 
adverse  pressure  gradient  flow;  note  procedure 
is  less  satisfactory  with  increasing  pressure 
gradient. 
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adverse  pressure  gradient  flows  respectively.  It  may  be  observed 
from  Figure  4.1  that  the  theoretical  profile  represents  the  con¬ 
stant  pressure  profile  well  in  the  overlap  zone  but  is  increasingly 
less  satisfactory  there  as  the  pressure  gradient  becomes  larger. 

It  may  be  observed  from  Figure  4.1  that  for  stations  8209  and 
8309,  a  larger  value  of  <  appears  to  be  indicated  in  order  to 
decrease  the  slope  of  the  profile  to  conform  to  the  data. 

However,  there  are  two  main  difficulties  associated'  with 
the  fitting  problem  and  it  is  worthwhile  to  discuss  these  here 
before  preceding  further.  The  first  problem  has  been  termed  a 
"low  Reynolds"  number  effect  by  Scharnhorst  (1978)  and  is  asso¬ 
ciated  with  a  failure  of  the  composite  analytical  profile  to 
adequately  delineate  various  regions  of  the  boundary  layer.  To 
understand  this  last  statement,  consider  profile  8309  in  Figure 
4.1;  the  solid  straight  line  is  apparently  tangent  to  the  curve 
in  the  overlap  zone  of  the  profile.  The  inverse  of  the  slope 
of  this  straight  line  may  be  obtained  graphically  and  is  indicated 
on  the  figure  as  k=0.266;  however,  the  value  of  k  used  to  obtain 
the  profile  in  the  figure  was  <=0.41.  The  reason  for  this  diffi¬ 
culty  may  be  clearly  observed  in  Figure  4.2  which  illustrates 
composite  similarity  profiles  for  two  values  of  the  Reynolds 
number  based  upon  the  displacement  thickness,  6*,  for  fixed  values 
of  6^  =  0,  K  =  0.016,  and  <  =  0.41.  In  preparing  this  figure. 


a  value  of  Re  5*  =  10^  was  chosen  since  it  is  typical  of  practical 
flows  that  occur  in  engineering  practice.  For  Re^*  =  10^  the 
range  of  y'^  that  exhibits  logarithmic  behavior  is  apparently 
relatively  short  and  the  slope  of  the  analytical  profile  in  this 
region  appears  to  be  greater  than  the  value  I/k.  The  reason  for 
this  behavior  is  associated  with  the  value  of  Re^*  which  relates 
the  inner  variable  y*"  to  the  outer  variable  n  by  y"*"  =  Re5*n, 

For  ROjj*  =  10^,  a  value  of  =  100  would  correspond  to  an  outer 
variable  value  of  n=0.1;  for  this  reason,  the  apparent  range  of 
logarithmic  behavior  in  Figure  4.2  for  Re^*  =  10^  is  relatively 
short  and  significant  portions  of  the  inner  and  outer  regions 
blend  together  in  the  overlap  zone.  To  show  the  effect  on  the 
profile  as  the  Reynolds  number  becomes  larger,  a  composite  simi¬ 
larity  profile  is  also  presented  in  Figure  4.2  for  a  value  of 
Reg*  arbitrarily  increased  to  10^,  It  may  be  observed  that  a 
logarithmic  region  emerges  over  a  wider  y^  range  and  that  the 
conposite  similarity  profile  corresponds  on  the  graph  more 
closely  with  the  input  log-law  behavior. 

A  second  difficulty  associated  with  the  fitting  of  the 
composite  profile  is  encountered  with  moderate  to  strong  pressure 
gradient  flows  where  8^  is  0(1).  It  may  be  observed  in  the 
series  solution  for  the  outer  profile  given  in  equation  (A. 4) 


i 


that  for  all  0,  the  series  contains  terms  of  the  form  a^-n^ 

log  n  in  addition  to  the  purely  logarithmic  term  (l/ic)logn;  two 
points  about  these  terms  are  germane.  In  the  first  place,  it 
may  easily  be  verified  that  such  terms  arise  in  the  series  because 
of  the  pressure  gradient  term  in  the  outer  layer  equation  and 
are  not  associated  with  the  choice  of  turbulence  model  per  se. 
Secondly,  although  logn  -»•  0  as  n  -*■  0  for  all  i  >0,  such  terms 
do  give  a  significant  contribution  for  n  /  0  and  the  purely 
logarithmic  behavior  of  (l/ic)log  n  will  only  be  realized  for 
very  small  n.  As  previously  indicated,  for  low  Re g*, very  small 
values  of  n  will  correspond  to  moderate  (but  not  large)  values 
of  y'*’  and  the  logarithmic  portion  of  the  analytical  profile 
becomes  obscure.  As  0.  increases  the  difficulty  becomes  more 
severe  since  it  may  easily  be  verified  from  equation  (A. 4) 
that  the  a^  become  increasingly  larger  as  increases.  This 
effect  further  causes  the  composite  profile  to  deviate  from 
logarithmic  behavior  for  small  n. 

The  failure  of  the  composite  profile  to  reflect  the  true 
input  logarithmic  behavior  for  low  Reynolds  number  on  a  graph 
is  of  some  concern  since  the  values  of  Reg*  which  are  character- 
istic  of  many  experiments  are  0(10  )  or  0(10  ).  To  attempt  to 
overcome  these  problems,  Schamhorst  (1978)  suggested  two 
approaches.  In  the  first  of  these,  a  full  three  parameter 
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optimization  of  S,  k,  and  K  was  carried  out.  Generally  this 
procedure  produces  very  good  representations  of  the  data; 
however,  as  Schamhorst  (1978)  points  out,  the  parameters  (S,k,K) 
tend  to  lose  physical  significance  in  such  method  in  two  ways. 

F’rst  the  input  value  of  <  used  to  produce  the  profile  is 
not  the  apparent  value  that  would  be  calculated  from  a  graph  of 
the  profile;  this  point  has  been  illustrated  in  Figure  4.1. 
Secondly,  Schamhorst  (1978),  upon  plotting  the  values  of  S 
obtained  for  profiles  with  various  values  of  observed  a 
trend  for  S  to  remain  approximately  constant  or  to  actually 
slightly  decrease  with  increasing  This  trend  is  opposite 
to  the  experiments  (Kline,  et  al.  1967)  which  suggest  that  S 
should  increase  with  increasing  b  . 

v» 

The  second  approach  was  suggested  by  Schamhorst  (1978) 
as  one  possible  way  of  overcoming  the  problems  of  the  three  para¬ 
meter  optimization  and  was  attempted  here.  The  main  ideas  are 
that  the  parameter  K  primarily  influences  the  quality  of  the 
fit  in  the  outer  region  while  the  parameter  S  mainly  influences 
the  inner  region;  on  the  other  hand,  influences  the  slope  of 
the  profile  in  the  overlap  zone  and  consequently  has  an  important 
effect  in  both  regions.  The  difficulties  associated  with  the 
strong  influence  of  varying  <  and  also  with  fitting  this  type  of 
profile  have  been  discussed  by  Schamhorst  (1978)  and  Weigand 
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(1978);  one  approach  found  to  be  acceptable  for  heat  transfer 
profiles  (Wiegand,  1978)  is  to  perform  a  one  parameter  optimiza¬ 
tion  on  K  holding  k  fixed;  once  the  optimization  had  been  com¬ 
plete,  a  new  value  of  <  is  obtained  from  a  graph  of  the  data  in 
inner  region  coordinates;  note  that  in  the  optimization  u^  will 
vary  and  hence  the  graph  of  u'*’  versus  y'*'  will  change.  This  graph¬ 
ical  iteration  process  starts  with  a  value  <=0.41  and  continually 
obtains  new  values  of  <  from  the  slope  of  the  data  points  in  the 
logarithmic  region  of  the  velocity  profile.  This  approach  is 
denoted  by  method  1  in  this  study  and  the  least  square  curve  fit 
error  results  of  the  first  iteration  are  presented  in  Table  4.1, 
First  iteration  curve  fits  for  the  last  data  station  of  the  zero 
and  favorable  pressure  gradients  of  Andersen,  et  al.  (1972)  are 
shown  in  Figure  4.1.  Unfortunately  method  1  failed  for  successive 
graphical  iterations  since  the  new  values  of  <  from  the  slope  of 
the  logarithmic  region  do  not  provide  improved  curve  fits;  in 
fact,  the  values  of  <  taken  from  the  data  give  a  lower  value  of 
K  than  0.41  although  a  larger  value  of  <  is  suggested  upon  com¬ 
parison  of  the  data  and  the  analytic  profile.  The  lines  drawn 
tangent  to  the  data  in  the  logarithmic  zone  are  depicted  in  figure 
4.1  as  broken  lines.  The  failure  of  this  procedure  is  again  due 
to  the  low  Reynolds  number  effect  for  large  In  the  next  sec¬ 
tion,  several  other  methods  are  discussed  in  an  effort  to  obtain 
improved  composite  profile  data  comparisons. 


4.2  Composite  Profile  Data  Comparisons 

Composite  similarity  profile  representations  were  obtained 
using  several  methods  in  which  certain  model  parameters  were 
optimized  using  the  direct  search  technique  while  holding  others 
constant.  The  five  methods  considered  in  this  study  for  profile 
representations  are: 

1.  Method  1  -  one  parameter  optimization  on  K  while  hold¬ 

ing  S  =  10.4965  and  k  =  0.41  constant; 

2.  Method  2  -  two  parameter  optimization  on  k  and  K  while 

holding  S  =  10.4965  constant; 

3.  Method  3  -  two  parameter  optimization  on  S  and  K  while 

holding  k  -  0.41  constant; 

4.  Method  4  -  three  parameter  optimization  on  S,  k,  and 

K; 

5.  Method  5  -  two  parameter  optimization  on  S  and  K  while 

holding  k  *  0,46. 

To  determine  which  method  offers  the  best  profile  representations, 
each  method  was  run  with  the  non -transpired  zero,  mild,  and  strong 
adverse  pressure  gradient  data  of  Andersen,  et  al.  (1972) (labeled 
8100,  8200,  and  8300,  respectively)  and  the  favorable  pressure 
gradient  data  of  Herring  and  Norbury  (labeled  2700  after  Coles 
&  Hirst,  1969).  These  data  sets  were  chosen  for  several  reasons. 
First,  in  all  data  sets  there  is  a  relatively  large  number  of 
points  in  the  inner  layer  and  this  is  important  in  assessing  the 
performance  of  the  unsteady  wall  layer  model.  In  the  second 
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place,  these  data  sets  reflect  a  wide  variety  of  pressure 
gradients  for  equilibrium  mainstream  velocity  distributions. 
Thirdly,  these  data  sets  are  relatively  recent  and  are  believed 
to  be  very  reliable.  The  least  squares  error  e  is  presented 
in  Table  4.1  for  each  of  the  five  optimization  methods.  The 
values  of  I  in  parentheses  represent  recalculated  values  which 
neglect  the  third  and  fourth  data  stations  for  the  Herring  & 
Norbury  (Coles  &  Hirst,  1969)  2700  flow  and  the  first  two  data 
stations  for  the  Andersen,  et  al.  (1972)  8100,  8200,  and  8300 
flows.  These  neglected  data  stations  generally  correspond  to 
upstream  stations  which  are  suspected  of  not  being  representative 
of  equilibrium  behavior.  The  recalculated  I  values  are  observed 
to  have  a  value  close  to  that  for  the  average  t  representing  all 
data  stations  except  for  the  2700  series  in  which  a  substantial 
improvement  is  noted. 

By  observing  the  curve  fit  error  for  the  mild  adverse 
pressure  gradient  (8200  series)  in  Table  4,1,  it  is  apparent 
that  substantial  improvanents  in  the  quality  of  the  fits  can 
be  realized  by  varying  more  than  one  parameter.  The  best  curve 
fits  are  obtained  with  methods  2  and  4  in  the  sense  that  the 
lowest  least  squares  errors  are  obtained  with  these  methods.  A 
similar  trend  may  be  observed  for  the  strong  adverse  pressure 
gradient  (8300  series).  Note  that  the  RMS  error  increases  in 
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the  range  50-70%  for  methods  1,  3  and  5  in  going  from  the  mild 
to  the  strong  adverse  pressure  gradient  case;  for  methods  2  and 
4,  the  percentage  increase  is  19  and  14  percent  respectively. 

The  8300  series  is  the  most  difficult  case  for  profile  represen¬ 
tation  because  of  the  low  Reynolds  nuirt)er  effects  and  pressure 
gradient  effects  which  have  increasing  importance  for  larger 
Sj.  as  discussed  in  §4.1.  On  the  basis  of  the  results  in 
Tables  4.1,  methods  1,  3,  and  5  may  be  ruled  out  as  effective 
parameter  adjusting  methods. 

To  further  assess  v4iich  of  the  surviving  methods  offers 
the  best  profile  representation,  the  complete  profile  optimiza¬ 
tion  results  for  methods  2  and  4  are  given  in  Tables  4.2  through 
4.9  along  with  the  corresponding  profile  representations  in 
Figures  4.3  through  4,10.  Evidently,  the  S,  k,  and  K  optimiza¬ 
tion  technique  of  method  4  has  the  lowest  average  root-mean- 
square  error,  e;  however,  this  result  is  not  totally  unexpected 
since  the  adjustment  of  all  three  parameters  offers  more  flexi¬ 
bility  in  curve  fitting  for  both  the  inner  and  outer  regions. 
Method  2  (k,  K  optimization  with  constant  S  =  10.4965)  offers  a 
close  second  choice  in  the  selection  of  a  profile  representation 
method  which  suggests  that  the  variation  of  the  S  value  does 
not  drastically  change  the  quality  of  the  curve  fits. 

An  important  observation  which  can  be  made  from  the  results 
of  the  three-parameter  fits  of  method  4  is  that  the  composite 
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Figure  4.3  Velocity  profile  comparisons  for  two  para¬ 
meter  ootimi ration  on  k  and  K  (with  S  » 
10.4965)  for  favorable  pressure  gradient  data 
of  Herring  and  Norbury  (1969^ 
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TABLE  4.4  Results  of  Method  2  two  parameter  fits  of  the  composite  similarity 
profile  to  the  mild  adverse  pressure  gradient  data  of  Andersen, 
et  al.  (1972). 


'1  LO  100  1000  lOOOO 

Figure  4.5  Velocity  profile  comparisons  for  two  parameter 
optimization  on  k  and  K  (with  S  ■  10.4965}  for 
mild  adverse  pressure  gradient  of  Andersen  et  al. 
(1972). 
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TABLE  4.5  Results  of  Method  2  two  parameter  fits  of  the  composite  similarity 
profile  to  the  strong  adverse  pressure  gradient  data  of  Andersen, 
et  al.  (1972). 


Figure  4.6  Velocity  profile  comparisons  for  two  parameter 
optimization  on  <  and  K  (with  S  =  10.4965)  for 
strong  adverse  pressure  gradient  of  Andersen 
et  al.  (1972). 
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Norbury  (Coles  and  Hirst,  1969). 
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TABLE  4.7  Results  of  method  4  three  parameter  fits  of  the  composite  similarity 
profile  to  the  constant  pressure  data  of  Andersen,  et  al.  (1972). 


Figure  4.8  Velocity  profile  comparisons  for  three  parameter 
optimization  on  S,  k  and  K  for  constant  pressure 
data  of  Andersen  et  al.  (1972). 
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Figure  4.9  Velocity  profile  comparisons  for  three  parameter 
optimization  on  S,  k  and  K  for  mild  adverse 
pressure  gradient  data  of  Andersen  et  al.  (1972). 
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TABLE  4.9  Results  of  method  4  three  parameter  fits  of  the  composite  similar 
profile  to  the  strong  adverse  pressure  gradient  data  of  Andersen, 


similarity  profile  encounters  difficulty  in  fitting  velocity  pro¬ 
files  measured  in  strong  adverse  pressure  gradient  flows.  An 
increase  in  e  with  pressure  gradient  can  be  observed  in  Table 
4.1;  this  is  not  entirely  unexpected  for  two  reasons.  First, 
true  self-similar  behavior  is  not  anticipated  in  situations  with 
large  adverse  pressure  gradients  since  boundary  layer  separation 
may  occur  before  similarity  is  achieved;  consequently,  it  may 
be  that  the  model  is  not  an  appropriate  one  for  large  adverse 
pressure  gradients.  In  the  second  place,  there  is  a  significant 
departure  of  the  composite  profile  from  logarithmic  behavior  in 
the  overlap  region  for  6^  0(1)  due  to  the  low  Reynolds  number 
effect  and  the  0^  effect  in  the  series  solution  as  discussed 
in  section  4.1.  The  effects  of  increasing  6  in  the  series 
solution  may  be  offset  by  increasing  x,  thereby  reducing  the 
magnitude  of  the  a^  and  b^.  coefficients  in  equation  (A. 4)  (see 
Appendix  A).  This  cancelation  effect  may  be  observed  in  the 
results  of  the  three-parameter  fits  in  Tables  4.6  to  4.9  which 
reveal  that  the  optimized  values  of  <  increase  with  As 
a  result,  it  may  be  concluded  that  the  three  parameter  optimiza¬ 
tion  process  adjusts  the  parameters  of  the  composite  similarity 
profile  to  correct  an  undesirable  low  Reynolds  number  effect 
or  large  6-  effect  in  the  series  solution.  However,  the  optimized 

V 

value  of  K  is  no  longer  associated  with  the  apparent  inverse 
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slope  in  the  logarithmic  portion  of  the  profile  as  taken  from 
a  graph  of  the  profile. 

To  summarize  the  investigation  of  all  the  optimization 
methods,  it  was  found  that  a  significant  reduction  in  error  of 
the  profile  fits  is  realized  by  optimizing  all  three  parameters 
of  the  profile,  namely  S,  tc,  and  K.  However,  the  three  para¬ 
meter  optimization  masks  any  physical  interpretation  originally 
associated  with  the  parameters  to  counteract  the  low  Reynolds 
number  and  large  effects.  Thus,  the  parameters  S,  ic,  and  K 
become  simply  profile  parameters  which  are  adjusted  to  obtain 
a  close  representation  of  measured  data.  The  objective  now  is 
to  obtain  correlations  for  these  parameters  based  upon  the  four 
equilibrium  flows  of  2700,  8100,  8200,  and  8300  in  order  to 
provide  the  basis  for  a  prediction  scheme. 

4.3  Parameter  Correlations 

In  order  to  develop  a  profile  prediction  method,  optimized 
parameter  values  were  obtained  for  the  three  parameter  optimiza¬ 
tion  of  method  4  and  the  two  parameter  optimization  of  method  2. 
Both  optimization  methods  used  the  negative  pressure  gradient 
flow  of  Herring  and  Norbury  labeled  2700  and  the  three  Andersen, 
et  al.  (1972}  flows,  8100,  8200,  and  8300  which  cover  a  range 
of  8^.  from  about  -0.4  to  1.5.  The  results  of  the  optimized 


data  fits  for  the  three  parameter  optimization  of  method  4  were 
then  plotted  versus  the  pressure  gradient  6.  as  shown  in  Figure 
4.11.  The  first  two  data  stations  of  flows  8100,  8200,  and  8300, 
and  the  third  and  fourth  data  stations  of  flow  2700  were  not 
included  in  the  correlations  because  of  the  observed  form  of  the 
experimental  velocity  profiles  which  did  not  appear  to  exhibit 
self-similar  behavior;  this  is  reflected  in  the  fact  that  the 
results  for  these  eight  data  stations  showed  a  substantial  devia¬ 
tion  from  the  straight  line  correlations  in  figure  4.11.  Prospec¬ 
tive  correlations  which  could  be  used  in  a  prediction  procedure 
were  obtained  by  fitting  a  least  squares  straight  line  and  quadra¬ 
tic  curve  through  the  optimized  parameter  values  represented  by 
the  solid  symbols.  The  resulting  correlations  for  the  three  para¬ 
meter  optimization  of  method  4  are  given  in  Table  4.10.  The 
RMS  curve  fit  error  for  the  correlations  show  that  the  quadratic 
curve  fits  give  only  a  very  slight  improvement  over  the  linear 
curve  fits.  The  three  parameter  optimization  correlations  indi¬ 
cate  that  K  is  nearly  constant;  most  of  the  variation  occurs  in 
K  which  in  effect  offsets  the  influence  of  in  the  series 
solution  for  the  outer  layer  mode.  A  slight  variation  occurs 
in  the  S  correlation  which  is  contradictory  to  the  experimentally 
observed  trend  for  the  dimensionless  time  period  between  bursts. 
The  correlation  indicates  that  S  increases  slightly  with  pressure 
gradient  while  the  visual  observations  of  Kline,  et  al.  (1967) 


6, 


FIGURE  4-n 


Parameter  correlations  obtained  from  three-parameter 
fits  of  the  composite  similarity  profile  to  four  sets 
of  data. 
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The  numbers  in  parenthesis  are  RMS  errors. 


show  that  S  decreases  with  increasing  pressure  gradient.  This 
contradiction  is  explained  by  recalling  that  the  three 
parameter  optimization  tends  to  mask  any  physical  significance 
that  may  originally  have  been  attributed  to  the  model  parameters. 
Another  point  which  concerns  the  value  of  S  is  worthy  of  men¬ 
tion;  in  the  wall-layer  the  data  exhibits  a  rapid  profile  vari¬ 
ation  which  may  give  rise  to  an  unnatural  weight  in  the  least- 
squares  error  calculation.  In  addition,  wall-layer  data  points 
also  have  been  shown  to  exhibit  velocity  measurements  which  are 
dependent  upon  the  pitot  probe  tip  geometry  (Andersen,  Kays 
&  Moffat,  1972).  All  of  these  factors  contribute  to  the  conclu¬ 
sion  that  the  slight  trend  in  S  with  pressure  gradient  obtained 
with  the  three  parameter  optimizations  may  not  be  significant. 

The  correlation  for  k  indicates  that  there  is  a  substantial 
influence  of  pressure  gradient  on  the  value  of  x;  this  is  partly 
due  to  the  influence  of  in  the  series  solution  which  is  offset 
in  the  optimization  by  increasing  the  k  value  for  larger 
Because  of  this  effect  and  the  Reynolds  number  effect  discussed 
previously,  the  values  for  k  obtained  from  the  correlations 
presented  in  Table  4. IQ  cannot  be  directly  associated  with  the 
inverse  slope  of  the  logarithmic  region  on  a  graph  for  profiles 
with  displacement  thickness  Reynolds  numbers  0(10  ).  The  curve 
fit  data  for  the  k  correlation  of  the  three  parameter  optimization 


exhibits  the  least  scatter  (lowest  RMS  error)  of  the  S,  k, 
and  K  correlations  which  tend  to  strengthen  its  reliability  for 
use  in  a  profile  prediction  scheme. 


As  a  point  of  interest,  correlations  based  on  the  two  para¬ 
meter  optimizations  are  also  included  in  Table  4.10.  It  may 
be  observed  that  these  correlations  show  the  same  trends  as  the 
results  based  on  the  three  parameter  optimizations  of  method  4. 
Note  that  the  dependence  of  K  on  6^.  is  weak  while  the  dependence 
of  K  on  6-  is  strong;  for  method  2  the  RMS  is  larger  than  for 
the  results  based  on  method  4.  Since  a  lower  RMS  is  also  observed 
in  the  actual  profile  fits  using  method  4,  this  procedure  is  con¬ 
sidered  somewhat  superior  to  method  2  and  the  correlations  in 
Table  4.10  based  on  method  4  are  recommended. 
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5.  MAINSTREAM  TURBULENCE 

5.1  Introduction 

The  influence  of  mainstream  turbulence  on  a  fully  developed 
turbulent  boundary  layer  flow  with  zero  pressure  gradient  has 
been  studied  by  a  number  of  investigators.  The  experimental 
investigations  of  Kline  et  al.  (1960),  Kestin  (1966),  Huffman 
et  al.  (1972),  Charnay  et  al.  (1971)  and  others,  indicate  that 
mainstream  turbulence  affects  the  turbulent  boundary  layer  velo¬ 
city  profile  in  many  ways.  In  particular,  a  thickening  of  the 
boundary  layer  with  increasing  mainstream  turbulence  level  is 
observed  along  with  a  progressive  increase  in  skin  friction. 

A  general  change  in  the  shape  of  the  non-dimensional i zed  mean 
velocity  profile  has  also  been  documented  in  which  there  is  a 
marked  reduction  in  the  "wake"  component  of  the  outer  layer  as 
the  mainstream  turbulence  level  increases.  Finally,  in  a  ther¬ 
mal  boundary  layer,  the  heat  transfer  at  the  wall  increases 
progressively  with  increasing  mainstream  turbulence.  To  con¬ 
sider  the  possibility  of  including  the  effects  of  mainstream 
turbulence  level  into  a  boundary  layer  prediction  method,  it 
is  appropriate  to  first  examine  the  turbulent  boundary  layer 
momentum  equation. 

The  continuity  and  momentum  equations  governing  nominally 
steady  turbulent  boundary  layers  have  been  given  in  equations 
(2.1)  and  (2.2)  in  connection  with  the  conventional  type  of 
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turbulent  boundary  layer;  these  equations  still  apply  when  the 
mainstream  is  turbulent.  Moreover,  the  boundary  conditions  in 
equations  (2.7)  and  (2.8)  for  the  mean  profile  are  still  correct. 
The  single  turbulence  input  to  the  equations  is  the  Reynolds 
stress  term  -u'v' ;  unfortunately  it  is  not  possible  to  incorpor¬ 
ate  the  fact  that  mainstream  turbulence  is  present  in  the  Rey¬ 
nolds  stress  other  than  through  a  correlation.  To  understand  the 
reason  for  this,  consider  the  turbulence  kinetic  energy  q^ 
defined  as. 


q2  =  +  v'^  +  (5.1) 

where  u'^,  v'^  and  w‘^  are  the  turbulence  intensities.  For  a 
laminar  mainstream  flow  ^  approaches  zero  at  the  boundary 
layer  edge  but  for  mainstream  turbulence 

^  ^  as  y  C5.2) 

The  mainstream  turbulence  level  Tu  is  defined  by 

where  Ug  is  the  local  mean  mainstream  velocity.  For  simplicity, 
assume  that  the  mainstream  turbulence  is  isotropic,  viz. 

u'^  =  v^  =  as  y  -*■  »  ;  (5.4) 

consequently. 
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v^,  -V  ...  ,  as  y  -t-  «.  (5.5) 

However,  u'v'  0  at  the  boundary  layer  edge  just  as  for  the 
normal  type  boundary  layer. 

It  should  be  remarked  that  the  momentum  equation  (2.2)  does 
not  contain  the  intensities  u'^,  v'^  since  these  terms  are  of 
lower  order;  even  if  these  terms  were  retained  in  equation  (2.2), 
the  dilenma  is  still  not  resolved  since  the  intensities  appear 
as  the  difference  (u'^-v'^)  which  still  must  vanish  at  the  boun¬ 
dary  layer  edge  for  isotropic  mainstream  turbulence. 

Since  it  is  not  possibl-j  to  incorporate  the  fact  that  main¬ 
stream  turbulence  is  present  in  the  boundary  conditions  for 
either  the  mean  profile  or  the  turbulence  terms,  the  other  pos¬ 
sibility  of  developing  a  correlation  for  the  eddy  viscosity 
model  parameters  is  investigated  here.  In  particular,  the  eddy 
viscosity  formula  in  equation  (2.55)  contains  the  parameters  K 
and  K  and  here  possible  trends  for  these  parameters  with  will 
be  considered.  In  addition,  a  possible  trend  in  the  inner 
region  profile  parameter  S  will  be  investigated.  Again,  this 
is  carried  out  by  a  systematic  adjustment  of  these  parameters 
to  obtain  a  best  fit  to  data  in  an  equilibriun  flow  but  now 
with  various  levels  of  mainstream  turbulence.  The  ultimate 
objective  here  is  to  provide  correlations  for  these  parameters 
which  could  be  used  in  a  prediction  method. 
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The  data  used  for  composite  profile  representations  was 
obtained  from  the  United  Technologies  Research  Center  boundary 
layer  wind  tunnel  (Blair  and  Werle,  1980).  Several  turbulent 
boundary  layer  mean  velocity  profile  data  sets  were  available 
for  zero  pressure  gradient  flows  with  mainstream  turbulence 
intensities  ranging  from  0.2%  to  6.4%.  These  data  sets  have 
a  large  number  of  points  that  provide  a  good  data  base  for 
profile  fits  and  represent  the  best  available  test  data  at  this 
time  which  reflects  the  effect  of  mainstream  turbulence. 

5.2  Composite  Profile  Data  Comparisons 

Composite  similarity  profile  representations  were  obtained 
for  two  methods  in  which  one  or  two  model  parameters  were 
optimized  using  the  direct  search  technique  while  holding  others 
constant.  The  first  profile  optimization  method  used  to  repre¬ 
sent  turbulent  boundary  layer  velocity  profile  data  with  main¬ 
stream  turbulence  was  a  one  parameter  optimization  in  K  with 
constant  values  S  =  11.025  and  k  -  0.44789.  These  constant  S 
and  K  values  correspond  to  the  values  taken  from  the  three 
parameter  optimization  correlations  in  §4.3  for  zero  pressure 
gradient.  Results  of  this  method  are  presented  in  Tables  5.1 
through  5.3  along  with  the  corresponding  profile  representations 
in  Figures  5,1  through  5,3;  note  that  the  18  data  stations  used 
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FIGURE  5.2  Velocity  profile  comparisons  for  one  parameter 
optimization  on  K  (with  S»11.025  and  ic*0. 44789) 
for  mainstream  turbulence  flow. 
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TABLE  5.3  Results  of  one  parameter  fits  of  the  composite  similarity  profile  to  mainstream 
turbulence  data. 


e  parameter 
d  ic»0. 44789) 


In  these  profile  comparisons  have  been  taken  from  the  study  of 
Blair  &  Wer1e  (1980) and  have  been  arranged  here  as  a  1000  series 
In  increasing  order  corresponding  to  Increasing  levels  of  local 
mainstream  turbulence.  Consequently  the  data  In  this  sequence 
represents  a  mixture  of  data  stations  from  the  four  basic  runs 
with  different  mainstream  turbulence  generators  which  were 
carried  out  by  Blair  &  Werle  (1980),  By  observing  the  least 
squares  error  e  tabulated  In  Tables  5.1  through  5.3  and  the 
velocity  profile  representations  presented  in  Figures  5.1 
through  5.3,  It  Is  apparent  that  the  quality  of  curve  fits  are 
acceptable  and  that  the  profile  model  can  be  used  to  demonstrate 
the  effects  of  mainstream  turbulence  on  turbulent  velocity  pro¬ 
files  with  good  success.  The  optimization  results  Indicate  that 
there  Is  a  trend  in  which  the  value  of  the  K  parameter  Increases 
with  Increasing  turbulence  level  Tu.  This  trend  was  anticipated 
In  the  sense  that  it  was  known  that  Increasingly  larger  main¬ 
stream  turbulence  levels  result  In  a  progressively  thicker  boun¬ 
dary  layer;  mathematically  larger  values  of  K  in  the  eddy  vis¬ 
cosity  formula  give  rise  to  a  thicker  boundary  Ijkyer, 

A  second  profile  optimization  method  was  considered  In 
an  attempt  to  obtain  Improved  velocity  profile  fits.  This 
method  was  Initiated  to  Investigate  possible  changes  In  both 
the  inner  and  outer  layers  of  the  velocity  profile  and  utilized 
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a  two  parameter  optimization  on  S  and  K  with  a  constant  value 
of  K  =  0.44789;  again  this  value  of  k  corresponds  to  the  value 
taken  from  the  three  parameter  optimization  in  §4,3.  The 
results  of  this  two  parameter  optimization  are  presented  in 
Tables  5,4  through  5.6  along  with  the  corresponding  profile 
representations  in  Figures  5.4  through  5.6. 

An  observation  which  can  be  made  concerning  the  results 
of  the  two  parameter  optimization  is  that  the  least  squares 
error  is  generally  smaller  than  for  the  one  parameter  optimiza¬ 
tion;  but  that  the  improvement  in  the  quality  of  the  curve  fits 
is  not  dramatic.  The  reason  for  this  behavior  may  be  explained 
by  close  examination  of  the  nature  of  the  curve  fits  near  the 
wall;  in  this  region  there  is  a  greater  curve  fit  error  for  the 
one  parameter  K  optimization  than  for  the  error  associated  with 
the  two  parameter  S  and  K  optimization.  The  remaining  sections 
of  the  model  velocity  profiles  for  the  one  and  two  parameter 
optimizations  are  similar  and  exhibit  almost  the  same  curve  fit 
error.  In  the  two  parameter  fit  the  value  of  S  adjusts  to 
minimize  this  error  for  the  data  points  nearest  the  wall. 
Unfortunately  these  data  points  are  usually  the  most  uncertain; 
it  is  also  important  to  note  that  the  S  parameter  values 
obtained  from  the  two  parameter  optimization  exhibit  a  great 
deal  of  scatter  and  any  attempt  to  derive  a  trend  from  this 
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FIGURE  5.4  Velocity  profile  comparisons  for  tw  parameter 
optimization  on  S  and  K  (with  k*0. 44789)  for 
mainstream  turbulence  flow. 
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information  was  considered  to  be  potentially  misleading.  Finally, 
the  magnitude  of  the  S  parameter  for  the  low  turbulence  inten¬ 
sity  cases  did  not  match  the  value  of  S  =  11.025  which  was 
obtained  from  the  S  correlation  of  §4.3.  All  of  the  above 
observations  suggest  that  the  one  parameter  K  optimization 
with  constant  values  S  =  11.025  and  k  =  0.44789  provides  the 
most  realistic  method  for  developing  velocity  profile  predic¬ 
tions.  The  objective  now  is  to  obtain  a  correlation  for  the 
K  parameter  based  upon  the  results  of  this  section. 

5.3  Parameter  Correlations 

In  order  to  develop  a  profile  prediction  method  optimized 
parameter  values  were  obtained  from  the  one  parameter  optimiza¬ 
tion  of  K  with  constant  values  of  S  *  11.025  and  ic  »  0.44789. 

This  optimization  method  examined  eighteen  zero  pressure  gradient 
data  sets  with  mainstream  turbulence  levels  which  cover  a  range 
of  Tu  from  0.002  to  0,0640.  The  results  of  the  optimized  data 
fits  for  the  K  parameter  were  then  plotted  versus  the  turbu¬ 
lence  level  Tu  as  shown  in  Figure  5.7,  Prospective  correla¬ 
tions  which  could  be  used  in  a  prediction  procedure  were 
obtained  by  fitting  a  least  squares  straight  line  and  quadratic 
curve  through  the  optimized  parameter  values  represented  by  the 
symbols.  The  resulting  correlations  for  the  one  parameter 
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K  parameter  correlation  obtained  from  one  parameter  fits  of  the  comp^lte 
similarity  profile  to  constant  pressure  test  data  with  mainstream  turbulence. 


optimization  on  K  are 


K  =  0.01489  +  0.39285  Tu  (.0008519)  (5.10) 

and 

K  *  0.01567  +  0.30941  Tu  +  1.35909  Tu^  (.0008462)  (5.11) 

for  the  linear  and  quadratic  fits  respectively.  Both  correla¬ 
tions  represent  Increasing  functions  of  mainstream  turbulence 
intensity  as  expected  and  discussed  in  section  5.1.  The  RMS 
curve  fit  errors  for  the  correlations  are  given  in  parenthesis 
following  the  equations.  These  RMS  values  indicate  that  the 
quadratic  curve  fit  gives  only  a  very  slight  improvement  over 
the  linear  fit;  however,  one  feature  of  the  quadratic  correla¬ 
tion  is  that  the  value  of  its  Intercept  is  close  to  the  value 
0.01591  obtained  from  the  three  parameter  correlation  of  §4.3 
for  turbulent  boundary  layers  affected  by  pressure  gradient. 

This  indicates  that  the  prediction  methods  for  pressure  gradient 
effects  outlined  in  §4  and  the  method  for  mainstream  turbulence 
are  compatible.  This  compatibility  adds  to  the  credibility  of 
the  correlations  for  use  in  profile  prediction.  To  obtain  a 
correlation  that  incorporates  the  zero  pressure  gradient  inter¬ 
cept  value  of  K  ®  0.01591  from  §4,3,  the  quadratic  curve  fit  of 
equation  (5.11)  was  rerun  with  the  intercept  held  at  0.01591. 

The  resulting  correlation  for  the  one  parameter  optimization 


6.  SUfflARY  AND  CONCLUSIONS 

In  the  present  study,  a  profile  model  for  the  mean  velocity 
in  a  nominally  steady  two-dimensional  turbulent  boundary  layer 
has  been  developed.  To  obtain  this  profile,  the  leading  terms 
in  an  asymptotic  expansion  for  high  Reynolds  number  for  the  mean 
velocity  were  derived  for  both  the  inner  and  outer  layers  of  the 
turbulent  boundary  layer.  A  self-similar  behavior  in  the  velo¬ 
city  was  assumed.  In  the  outer  layer,  a  simple  eddy  viscosity 
formula  was  assumed  containing  two  parameters  ic  and  K;  the  outer 
layer  self-similar  profile  satisfies  an  ordinary  differential 
equation  which  was  solved  numerically.  For  the  inner  layer,  an 
analytical  profile  model  which  is  based  on  the  observed  coherent 
structure  of  the  time-dependent  flow  in  the  wall  layer  was  used; 
this  model  contains  the  parameters  k  and  S.  A  composite  profile 
which  spans  the  entire  boundary  layer  was  defined  and  this  pro¬ 
file  is  in  general  a  function  of  the  three  parameters  (K,  k,  S). 

A  computer  code  was  developed  for  which  any  or  all  of  the  three 
profile  parameters  may  he  varied  to  obtain  a  best  fit  to  a 
given  set  of  experimentally  measured  profile  data;  this  code 
is  described  in  Appendix  C,  A  test  case  for  the  code  is  pre¬ 
sented  in  Appendix  D, 

There  are  two  potential  uses  for  this  optimization  procedure. 
In  the  first  of  these,  a  very  close  representation  of  a  given 
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set  of  velocity  data  may  be  obtained  by  carrying  out  a  three 
parameter  optimization  using  the  code  described  in  Appendix  C 
even  in  situations  where  the  flow  does  not  exhibit  self  similar 
behavior.  The  second  use  is  rather  more  significant  and  is 
associated  with  the  development  of  the  basic  model  for  a 
boundary- layer  prediction  method.  By  carrying  out  a  series  of 
optimization  studies  for  given  sets  of  data,  it  is  demonstrated 
in  §4  and  §5  that  trends  in  the  optimized  values  of  Cs,k,K)  may 
be  observed;  from  chese  trends  correlations  for  the  model  para¬ 
meters  for  a  ph/‘i'cal  effect  may  be  obtained.  In  the  present 
report,  two  such  :t;dies  have  been  carried  out  and  correlations 
have  been  developed  for  the  effects  of  both  pressure  gradients 
and  mainst’"eam  t’jrbulaice. 
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APPENDIX  A 


OUTER  LAYER  SIMILARITY  SOLUTION 
A. 1  INTRODUCTION 

In  this  section,  the  outer  layer  similarity  equation 
A  dU,  dU, 

^  +  (H2e^)n  -3^  +  23^U^  =  0  .  (A.l) 

is  solved  to  obtain  the  outer  layer  velocity  profile.  Because 
of  the  irregular  logarithmic  behavior  near  ii=0,  two  types  of 
solution  methods  >  'e  used  and  these  consist  of  a  series  solu¬ 
tion  for  small  n  and  a  numerical  solution  for  large  n.  The 
matching  point  where  the  series  and  numerical  solutions  are 
required  to  merge  smoothly  is  taken  to  be  »  K/2ic, 

A. 2  SERIES  SOLUTION  FOR  SMALL  n 

The  eddy  viscosity  function  e(n)  was  selected  such  that 
eCn)  approaches  a  linear  variation  Kn  exponentially  quickly  as 
n  approaches  zero.  Consequently  for  small  n»  the  outer  layer 
similarity  equation  (A,l)  becomes 

d^U,  dU, 

SiF  *'  -ar  “1  '  “  • 

For  equation  (A, 2),  a  series  solution  of  the  form 


00 


u 


1 


i 

n  , 


{A.3) 


Is  assumed  and  the  indicial  equation  yields  a  double  root  a=0. 
The  solution  thus  takes  the  form  of  a  regular  power  series 
plus  a  logarithmic  term  of  the  form, 

U,  =  1  b„  n*'  +  {  2  a  n"’}  log  n  ,  (A. 4) 

‘  n=0  "  n=0  " 

with  derivatives 

00  00  00  n*7 

E  h  b  n”  +  {  E  n  a  n  }log  n  +  s  n 
n=l  ”  n=l  "  n=0  ’ 

{A.5) 

E  n{n-l)b_n''’^  +  {  r  n(n-l)a_n'''~^}log  n 
n=2  "  n=2  " 

+  E  n  a_n"”^  - -T  ^  a„{n-l)n"‘^  .  (A. 6) 

n=l  "  "I  n*2 

As  TT+O,  the  required  behavior  for  is  given  by  equt^ton 
(2.63)  and  is 

U-,  -  —  log  n  +  C  ;  (A. 7) 

I  1C  0 


dU^ 

dn 


and 


d^U^ 
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therefore 


a 


0 


1 

K 


and 


(A.8) 


Recursion  relations  for  a^^  and  are  obtained  by  substitution 
of  and  its  derivatives  from  equations  (A. 4),  (A. 5),  and  (A. 6) 
into  equation  (A. 2).  It  may  easily  be  shown  that 


-{2s  +n(l+2eJ}a„ 
^n+1  " - - 

and 


(A.9) 


+  (2Sj.  +  n(l+26j.))b^} 


(A.IO) 


A.  3  NUMERICAL  SOLUTION  OF  OUTER  REGION  FOR  LARGE  n 

A  numerical  solution  to  the  outer  layer  similarity  equation 
(A.l)  is  calculated  in  the  range  n^,  <  n  <  where  is  the 
matching  point  to  the  series  solution  and  represents  the 
outer  bound  of  the  numerical  mesh  which  is  chosen  large  enough 
to  ensure  no  significant  change  in  the  solution.  In  order  to  obtain  an 
accurate  numerical  solution,  a  small  mesh  size  is  needed  near 
the  wall  whereas  far  from  the  wall,  the  solution  decays 
rapidly  and  a  larger  mesh  may  be  used.  For  this  reason,  a 
numerical  method  developed  by  Kellor  (1969)  which  permits  a 
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where  -j  is  the  variable  mesh  length.  Equations  (A.  14) 
and  (A. 15)  are  approximated  at  the  midpoint  of  nj  and 
using  a  central  difference  for  the  first  derivative  and  aver¬ 
age  values  for  u  and  v.  As  a  result,  equations  (A. 14)  and 
(A.  15)  become* 


— =  —  (v-  *  +V") 
2  ^’'j+1 


and 


(A.17) 


j 

where 

’’iH  ■  ’JH  "  ‘’'"j  *  • 

The  auxiliary  variable  v^^^  is  eliminated  using  equation  (A.17) 
to  obtain 


A  similar  procedure  is  used  to  approximate  (A. 14)  and  (A. 15) 
at  the  midpoint  of  Tij_i  and  nj.  Eliminating  the  variable  Vj._^ 
in  a  similar  fashion  reveals 


3- 


(A.20) 


h“  (Uj-Uj.i)  = 


By  adding  (A. 19)  and  (A.20)  to  eliminate  v.,  it  may  be  shown 
that  the  finite  difference  equations  reduce  to  the  general 
form 

"j  “j+l  *  "j“4  *  '  ’*4  • 


where 


a  • 
J 


d,  =  0  , 


This  general  form  can  be  used  to  generate  a  syston  of  N-1 
finite  difference  equations  for  a  mesh  with  j=l,2,3,...,N 
number  of  grid  points.  The  system  of  equations  forms  a 
tri-diagonal  matrix  which  can  be  solved  in  a  direct  and 
efficient  manner  by  using  the  Thomas  algorithm;  this 
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algorithm  defines  two  functions  5^  and  F.  In  which  the  system 

J  J 

(A. 21)  may  be  solved  as 


^j+1 


(A.22) 


where 


'1* 


(A.23) 


for  j=l,2,3,...,N-l.  After  the  arrays  F.  and  5-  are  calculated, 

w  J 

the  solution  u^  is  obtained  from  (A.22)  by  back  substitution 
from  j=N-l  to  j=l. 


A. 4  SOLUTION  PROCEDURE 

The  procedure  used  in  this  study  to  obtain  a  smooth  match 
between  the  series  and  numerical  solutions  at  the  match  point 
is  outlined  below.  Let  the  solution  to  equation  (A.l)  be 
represented  by 

J  u  ,  0  <  n  <  rim 
u  =  <r 

I  “n  *  ^  1  <  ^0 

where  u^  and  u^  denote  the  series  and  numerical  solutions 
respectively.  First,  arbitrary  value  for  the  constant 
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.  From  this  value. 


in  equation  (A. 7)  is  guessed,  say 
a  series  solution  may  now  be  calculated  for  0  <  n  <  say 
Ug^^^.  In  particular,  the  series  solution  is  used  at  the 
match  point  to  obtain  and  its  derivative  du5^^^(nnj)/dn. 

Using  the  numerical  solution  for  ji  n  is 

calculated,  say  u^^^^(n);  the  derivative  du^  ^(njjj)/dn  is  cal¬ 
culated  numerically  using  a  six  point  forward  difference  for¬ 
mula  (Abramowitz  4  Stegum  1972,  p.  914).  It  is  worthwhile 
to  note  that  at  this  step  the  resulting  solution  will  not  have 

a  continuous  derivative  at  n=n^.  A  second  arbitrary  value 

(2l  (2)  (2) 

for  Cq  is  guessed,  say  Cq'  ‘  from  which  '  and  u^^'  '  are 

calculated  using  a  similar  procedure  as  that  stated  above. 

Again,  the  resulting  solution  will  not  have  a  continuous 

derivative  at  n*n^. 

The  final  solution  with  a  continuous  first  derivative  at 
n=n^  is  a  linear  combination  of  the  solutions  previously  found. 
Thus, 

Us  =  B^Us^^^  +  B2Us^^^  for  0  <  n  <  n^,  (A. 24) 

and 

“n  *  ^o''  %  1  "  i  .  (*-25) 

where  B-j  and  B2  are  constants.  One  condition  Imposed  on  B^ 
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and  B2  Is 


(A.26] 


+  B2  ®  1  , 

in  order  that  u  satisfies 

u  -  —  log  n  +  as  n  -►  0  . 
K  ^  0 


A  second  condition  relating  B-j  and  B2  is  that  the  first 
derivatives  of  equations  (A. 24)  and  (A. 25)  must  be  equal  at 
the  match  point  according  to 


du. 


(1) 


du. 


(1) 


du 


(2) 


dn 


dn 


+  B2{- 


du. 


(2) 


~S^ 


dn 


»  0. 


n®n 


m 


(A.27] 


Equations  (A.26)  and  (A.27)  can  be  solved  to  yield  the  values 
of  B-j  and  B2.  The  true  value  of  Cq  is  thus 

Co  =  .  (A.28] 


In  this  manner,  a  solution  to  equation  (A.l)  is  obtained  con¬ 
sisting  of  a  series  solutim  for  n  <  n,P  and  a  numerical  solu¬ 
tion  for  n  >  n^;  this  solution  and  its  first  derivative  are 
continuous  at  the  matching  point  n=n^.  Since  equation  (A.l) 
is  a  second  order  equation  all  derivatives  are  therefore  con¬ 
tinuous  at  n®%,. 
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APPENDIX  B 


THE  H  FUNCTION 


The  H  function  is  defined  as 


H  (n)  =  j 


e~^^dt  dy  dx  , 


(B.l) 


which  satisfies  the  differential  equation 


+  2n  H  =  ^  erf  (n) 


{B.2) 


The  following  expansion  is  readily  obtained 


H(n)  = 


_  e 


imv: 


T 


(B.3) 


where 


a(j)  =  a(j-l)  +  j  ,  a(l)  =  1  , 


(B.4) 


which  is  uniformly  convergent  for  all  n.  An  asymptotic  expan¬ 
sion  used  to  evaluate  H(n)  as  n  -*■  ®  is 

E(n)  -#{logn+^-i  ^  >  *  (B.5) 


where  is  Euler's  constant  equal  to  0.57721566...  . 
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APPENDIX  C 

THE  OPTIMIZATION  CODE 


oonooooooooooooooooooooooooon 


PROGRAM  HAIM (INPUT .OUTPUT, TAPE5aIM»UT«TAPE6»0UTPUTI 
COMNON/VAR/IVAROI  .XITtST 
EXTERNAL  F 

•  •  •  THIS  IS  THE  HAIM  PROGRAM  WHICH  INITIATES  THE  DATA  FITTING 

•  •  •  BY  CALLING  SUBROUTINE  PROFIT. 

•  •  •  NPAR  -  NUMBER  OF  PARAMETERS  DESIRED  FOR  THE  PROFILE  FIT. 

•  •  •  REFERENCEI  AN  OPTIMIZATION  TECHNIQUE  FOR  THE  DEVELOPMENT 

•  •  •  OF  A  TWO-DIMENSIONAL  TURBULENT  BOUNDARY  LAYER 

,  ,  ,  MODEL,  YUHAS,  L.J.  MASTER'S  THESIS 

•  •  •  LEHIGH  UNIVERSITY 

•  •  •  ITERH  -  MAXIMUM  NUMBER  OF  CYCLES  TO  BE  PERFORMED  IN  ORSRCH 

(NORMALLY  50) 

•  •  •  IPRINTsl  IF  PRINTING  OF  INTERMEDIATE  RESULTS  IN  ORSRCH  IS 

DESIRED  DURING  THE  COURSE  OF  THE  OPHNIZATION. 

>0  IF  NO  PRINTING  IN  ORSRCH  IS  DESIRED. 

•  •  •  EPS  -  THE  EXIT  TOLERANCE  IN  ORSRCH  IF,SAY,EPS«l.E-5  THE 

ITERATION  MILL  CONTINUE  IN  ORSRCH  UNTIL  THE  PARAMETERS 
BEING  OPTIMIZED  HAVE  CONVERGED  TO  5  SIGNIFICANT 
FIGURES 

THE  PARAMETERS  CAPK.S.KAPPA  TO  BE  OPTIMIZED  ARE  INITIALIZED 
IN  THE  INPUT 

IVARtDa  0  NOT  OPTIMIZED 
«  1  OPTIMIZED 
MHERE  I  s  1  CAPK 
=  Z  S 
S  3  KAPPA 

REAO(5,10) (IVAR(I) .Isl.S) 

NPARsQ 
DO  30  1=1,3 
30  NPAR=NPAR»IVAR(I) 

READ (5,201  TTERN, IPRINT.EPS.NA 

CALL  PROFIT(NPAR,ITERM,IPRINT,EPS,NA) 

10  F0RNAT(3I2) 

20  FORMAT  t2I(>,E6.1,I<»l 
STOP 
END 


ouooouuuuu  ouuu 


SUBROUTINC  PROFIT! Nf>AX«ITERN*IPIIINr,EP5*N«) 

OINENSION  ITITLEISl«LUNITI4l»XFC3l 

OIHENSION  LSYS<2I«I0ENI15}  tXdSI  «UE(19>  .UTAUIISI  tOUEDXIlS)  .PPLUSCl 
»5I «0ELTAtl61 fOELSTRt IS)* theta: 15) .SHAPE  CIS) «EPSIC15) .OELTCIS) 
OIHEHSION  PPt.USE(15)  .UTAUEC IS)  .ESI  IS)  .TOPLUSCIS)  *C(15)  tSETEClS) 
OIHENSION  CF(1S)*CFE(15) *CI ( IS) tCAPKC IS) .XKAPCIS) «8T CIS) . TUT  CIS) 
OIHENSION  FSC3)«DXC3)«0XERRC3) *0XFC3) ,OXERFC3) ,XNC IS) 

CO HNON/  OAT A/NP  Aft  »N OP  *  XNU  «  U I  *  OU lOX  ,  r 0  C  90  )  « UO  C 90  )  *  OL  STR* UA  C  so  )  « 
»YOPL(ia001 tUAPLClOOO) «NPPT 
C0HH0N/VAR/IVARC3) *X0  C3) 

COHHON/UOUT/NO.ETA 1^00) «  OETA  CfcOO) *UC  AOO ) »ETAH*  ETAO  «COUT«CIN« 

1 S  •  TNOTP  .  OE  C.  I N  «  OELOUT  •  UST  Aft  «  BET  AC 
EXTERNAt  F 

OATA  t.UNIT/2HFT*2H  N*2HIN*2HCH/tLSTS/7HEN6LISH ,7H  HETRIC/ tCBAR/0./ 
OATA  H6HT/0.1/ 

OATA  NPPT/400/ 

Na«NA 

•  •  •  SUBROUTINE  PROFIT  READS  IN  THE  EXPERINENTAL  OATA  AND 

•  *  •  INITIATES  THE  OPTIMIZATION  AT  EACH  OATA  STATION  BY 

•  •  •  CALLING  SUBROUTINE  ORSRCH 

•  •  •  INPUT  IDENTIFICATION  FOR  OATA  RUN 
ID  •  A  FOUR  DIGIT  NUMBER  CUSER  SUPPLIED)  TO  IDENTIFY  OATA 
lUNIT  >  ZERO  FOR  ENGLISH  UNITStONE  FOR  METRIC  UNITS 
ITITLE-  TITLE  OF  OATA  SET  CSO  CHARACTERS  HAXINUH) 

REAO(S*100)IOtIUNIT, (ITITLECI) tlslvS) 

IUP1»IUNIT*1 
IUP3»IUNrT»3 
CONV«12. 

iFCIUNIT.EQ.DCONValOO. 

NPARsNPAX 

•  •  •  XNU  -  KINEMATIC  VISCOSITY 

•  •  •  NSTA  -  NUMBER  OF  DATA  STATIONS 

REAOIStlOZ)  NSTA 
C 

C  INPUT  PARAMETERS  FOR  PLOTTING 
C 

C  •  •  •  NCYC  -  NO.  OF  X  CYCLES 

C  •  •  •  NPLTPP  -  NO.  OF  ZERO  PLOT  TICKS  ON  Y  AXIS 
C  •  »  •  XLC  -  LENGTH  OF  EACH  CYCLE  ON  X  AXIS 
C  •  *  •  INCPL  -  INCREMENT  LABELS  ON  Y  AXIS 
C  •  •  •  NL  -  NO.  OF  NON-ZERO  TICKS  ON  Y-AXIS 

C  •  •  •  XL  -  LENGTH  BETMEEN  TICKS  ON  Y-AXIS 

C 

REAOC5«90)  NCVC«NPLTPPtXLC«XL*NL*I9CPL 
90  FORNATC2I4*2F10.2«2II») 

00  10  Nsl.NSTA 

ID»IO*l 

IOENCN)-IO 

HRITEC6«103)  CITITLECD*  I«1«S)  *10 
C 

C  •  •  •  X  -  LOCAL  VALUE  OF  X-LOCATION  ON  THE  HALL 

C  •  •  •  UE  -  LOCAL  MAINSTREAM  VELOCITY 
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C  *  *  *  UTAUE  -  ESTIMATE  OR  EXPERIHENTALLT  QUOTED  VALUE  OF  UTAU 
C  •  •  •  DELTA  -  EXPERIMENTAL  VALUE  OF  BOUNDARY  LAYER  THICKNESS 

C  »  •  •  OELSTR-  EXPERI*NTAL  VALUE  OF  DISPLACEMENT  THICKNESS 

C  •  •  •  THETA  -  EXPERH€NTAL  VALUE  OF  MOMENTUM  THICKNESS 

C  •  •  •  DUEDX  -  LOCAL  VALUE  OF  THE  MAINSTREAM  VELOCITY  GRADIENT 

C 

REAO!5«IOi)XNU 
XN(Nl3XNU 

RE AD<5, 104 > X <N} «UEtN>«UrAUElN>, DELTA (N> .OELSTR INt , THETAiN) .OUEOXIN 
♦) tSETEINI 
READXS.130I  TUTCNI 

•  •  •  CFE  -  EXPERIMENTAL  VALUE  OF  THE  SKIN  FRICTION  COEFFICIENT 

•  •  •  SHAPE-SHAPE  FACTOR 

CFErNI=2.»  :UTAUErNl/UEtMII**2 
SH APE f N > 30ELSTR ( N1 /T HETA (N ) 

«  •  •  IF  AN  ESTIMATE  OF  UTAU  IS  NOT  AVAIL ABLE.THE  VALUES  OF 

•  •  *  UTAU  ARE  TO  BE  READ  IN  AS  ZERO  AND  AN  ESTIMATE  FOR 

•  •  «  UTAU  IS  COMPUTED  USING  THE  LUONEIG- TILLMAN  CORRELATION. 

•  •  •  NOTE  THAT  THIS  IS  ONLT  USED  AS  A  STARTING  ESTIMATE 

•  •  •  FOR  UTAU  IN  THE  OPTIMIZATION  PROCEDURE. 

IF(ABSIUTAUE(NI).GT.i.E-16tG0  TO  15 

CFE(NI*0.246*<OE(NI*TMeTAIN»/<XNU*50NV»»**<-0.268»*18.**<-.6r8*SHI 
♦PE<NII 

UTAUE (Nt SUE tNl«SQRTi0.5PCFEfNn 
15  CONTINUE 

•  •  *  PPLUSE-EXPERIHENTAL  value  OF  THE  INNER  REGION  PRESSURE 

•  •  •  GRADIENT  PARAMETER 

PPLUSE (N)s.XNU*UE( Nl sOUEOX  INI/ lUTAUE INI »*3) 

MRITE(6« 1051 X  INI .LUNITIIUPII .UEINI  .LUNITIIUPII 

NRITEIG.lOeiOUEDXtNI .UTAUEtNI .LUNITIIUPII .DELTAINI tLUNIT(IUP3l «OEI 
»STR(NI .LUNITIIUPII .THETA  INI .LUNITI IUP3I 

•  •  •  YD  -  Y  LOCATIONS  OF  THE  EXPERIMENTAL  DATA  POINTS 

•  •  •  UO  -  EXPERIMENTALLY  MEASURED  MEAN  VELOCITY  AT  YD 

SET  UP  AXIS  FOR  PLOTTING 

IFIN.GT.NPLTPPI  go  to  26 
IFIN.GT.ll  GO  TO  25 
CALL  PLOTtXLC*NCYC»2..2.0.-3l 
CALL  LOGAXINCYC.XLC.HGHTI 
CALL  SHFTYAXINL.NPLTPP.XL.INCPL.HGHTI 
CO  TO  22 
26  IF:N.GT.NPLTPP«^1I  go  to  25 

CALL  PLOTIXLC*MCYC«’2..-KL*IMPLTPP-ll,-3l 
CALL  LOGAXINCYC.XLC.HGHT* 

CALL  SMFTYAXINL.NSTA-NPLTPP.XL.INCPL.HGHTI 
GO  TO  27 

25  CALL  PLOT(O.O.XL.-3I 
27  CONTINUE 

REAO(5.102I  NOP 
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RCA0(5«99)  XOCll *Xa(2l,Xa(3) 

99  F0RMAT(3Fia.a) 

REAa;9«107l  CVO(I) *Isl«NOP) 

REAOIS.iar)  (UDttI  ttsl.NOP) 

MRITE(6«10ai  NOP *1. UNIT (IUP3) 

WRITE (6. 1091  (VO(I>.UO(I),I=l«NDPI 
□0  20  la  1*  NOP 
20  TO  (I)=yO<I)/CONW 

DELTA (Nl aOELTA (Nt/CONV 
aLSTR=OELSTR(NI/CONV 
USTARaUTAUE(Nl/UE(N> 

UIaUE(NI 
OUIOXaOUEOX<NI 

INITIALIZE  STEP  SIZES  FOR  ORSRCH 

OXdla.OOl 
OX(21a0.5 

ox:3)ao.oi 

OXERRdlaO. 00001 
OXERR(2ts0.0001 
OXERR(3)a0.0001 

ARRANGE  VARIABLES  BEING  OPTIMIZED 
IN  ORDERED  FORM  FOR  OPTIMIZATION  IN  ORSRCH. 


JCaO 

00  30  Iat,3 

IFdVARd)  .EQ.OtGO  TO  30 
JCaJC^l 
XF(JCIaXOdl 
OXFUCIaOX  dl 
OXERFUCMOXERRd) 

30  CONTINUE 

INITIATE  OPTIMIZATION  IN  ORSRCH 

CALL  0RSRCHfNPAR,3»FfFF,FS.XF,0XF,0XERF,ITERM,EPS» ITER. 
ilPRINT.IER) 

OPTIMIZATION  COMPLETE.  REARRANGE  VARIABLES  IN  CORRECT  ORDER. 


JC=l 

00  <»0  Ialf3 

IFdVARdI.EQ.  0)GO  TO  bO 
X0d)*XF(JCI 
OX  d)sOXF»JC> 
OXERR(I)aOXERF(JC> 
JCaJC»l 
40  CONTINUE 

EPSHNlaSQRT(FF» 

EBARaEBAR»rPSI(N) 

DELTA  IN) aOELTA(N)«CONV 
CAPKINIaXO  111 
ES  IN)aX0  12) 

XKAP(N)aX0<3) 

utau<n)*ustar*ue;n) 


122 


o  n  o  n  n 


CF(Nls2.*(UTAU(N)/ue(N)t **2 
CI(NI=CINF(EStNI ,TaPLUSt Nl «X0 :3) «0 . • 1. ) 

BT(N)3BETAC 

MRITE(6«liai  ITER.CPSKNIvNPAR 
MRITE(6tlll)  (Xa(JI«Jsl*31 vUTAUtNt 
HRITE16.112ICI:NI.T0PLUS!N) 

NRITEI6«1131  BETAC 

C  HATCH  CONDITION  CRITERIA  FOR  INNER  AND  OUTER  REGION  VELOCITY  PROFILES 
HRITE(6fll4l 

HRITE(6,115)  ETAH.NO.ETAO 
C  INNER  AND  OUTER  LENGTH  SCALES 
MRITE(6«2Q0) 

HRITE(6.201I  OELIN.OELOUT 
HRITEt6«116> 

00  50  1*1. NOP 

YP*Y0<I1/0ELIM 
UOPxUOdl/USTAR 
UAPsUACII/USTAR 
UOPHUAP>UOP-UAP 

50  WRITE  (6,117)  YP,UO(l),UOP,UAtI),UAP,UOPNUAP 
SET  UP  FOR  PLOTTING 

COHPUTE  ENOUGH  PLOT  POINTS  FOR  A  SMOOTH  CURVE 
FOR  ETA  LESS  THAN  ETAN 

NNPTS*A00 

ETPL*  ET  AM- 1 , "OEL IN/0  EL OUT 
OELPaETPL/ FLOAT (NNPTS) 

VOPL(1I»1.000«OELIN 
00  55  Is 1, NNPTS 
IFd.EQ.l)  GO  TO  54 
YOPL  d)  a YOPL  d-1 1  ♦OELP*OELOUT 

54  VOL=YOPLd)/OELIN 
ETAD=YOPLd)  /OELOUT 

CALL  U1SER(XKAP<NI  , BETAC, COUT, ETAD,U1 ,U1D, EPS) 

HAKEaUl-ALOG (ETA  0) /XKAPf  N) -COUT 

UAPL(I)3USTAR* (UP(TOL,S,TNOTP,CIN,XKAP(N),0.) »NAKE) 

55  CONTINUE 

00  56  1*1, NNPTS 

YOPL d)a (FLOAT (NCYCI*XLC)/ALOG(10.**FLOATINCYC))*(ALOG(YOPLtI» /DEL 
IIN)) 

UAPL d )*XL«NL«UAPL (I ) /(NL«INCPL*USTAR) 

56  CONTINUE 
C 

C  FOR  ETA  greater  THAN  ETA  NATCH,  USE  THE  NUMERICAL 
C  SOLUTION  ALREADY  CALCULATED 
NUAaNNPTS^NO 
NEXaNNPTStl 
00  55  I*NEX,NUA 
YOPLd)«eTAd-NNPTSI*OELOUT 
YOL*YOPL(I)/OELIN 
ETAOaETA  CI-NNPTS) 

NAKEaU(l-NNPTS)-ALOC(ETAO)/XKAP(N)-COUT 
UAPL (I )«USTAR« (UP(TOL,S, TNOTP.CIN, XKAP( N) , 0 .) »MAKE) 
YOPLd)«(FLOAT(NCYC)  «XLC) / ALOG (10. •»FL0AT (NCYC ) )  •  t  ALOGCYDPL  d) /OEL 
tIN)) 

UAPL(I)»XL*NL*UAPLd)/(NL*  IMCPL'USTAR) 


59  CONTINUE 

CALL  PL0T(Y0PL(i)*UAPLfl)*3> 

00  60  laZtNOP 

CALL  SYMBOL ( (FLOAT (NCYC) 'XLC) /ALOS:iO .••FLOAT ( NCYCI I •lALOG tYOt II /O 
1 EL IN) I «  XL*NL«UO{ I) /( NL*INCPL*USTAP) 

60  CONTINUE 

CALL  PLOT(Q.,0.«3) 

XLIM-(FLOAT(NCYC)*XLCI/AL06(10.««^LOAT(NCYCn*  :ALOGtYO:NOP)/OELIN) 

»» 

00  65  Isl.NUA 
NUAL^I 

IFCVOPLd)  .GT.XLIH)  GOTO  66 
66  CALL  PLOTtYOPLCt)«UAPL(I),2l 
66  CONTINUE 

XIOsYOPL(NUALI •2.*HGHT 
YI03UAPL(NUAL)-HCHT/2. 

CALL  NUHBEPtXIOf YID,HGHT«I0EN(NI«a.t2Ht4) 

10  CONTINUE 

EBAR«EBAR/FLOAT(NSTA) 

MRITE (6*110)  (inTLE(I>,lsl«5>,LSYS(ZUPl> 

00  70  1>1«NSTA 

70  MPITE(6*119) IOEN(I)»X(D •UE(I) .UTAUE(I) .XN ( I) « BETE (Z)*OELTA (X) * 

•  OELSTRdItTHETAd)  «TUTa) 

HRITE<6*120I  NPAR 

00  BO  Z«1«NSTA 

BO  MRITE(6«121)I0EN(It«8T»),UTAU(II*ES(XI  »CAPKdl*XKAPd>tEPSZri)* 

•  TOPLUSdltCKII 
WRITE (6*122)  NSTA*E8AR 

100  FORMAT  :2I4*SA10) 

101  FORMAT(E10.3) 

102  FORMAT (ID 

103  FORMAT(1H1«5X*5A10*X  >  F*I6) 

106  FORMAT:7F10.a*F10.2) 

105  F0RMAT(/*6X*«X  «  «*F6.3«1X*A2*6X*AUE3  A *F7. 2*1 X* A2 *P/SEC7 ) 

106  FORMAT l/*16X*AEXPERIMENTAL  VALUES* * 7 * 14X*A0UED A*F7.3**  1/SEC#*7 
l*16X,«UTAUa  A* F7 .3 *1X*A2 *A/SECA*/* 16X* AOELTA«  A*F7 .6*1X* A2 */*l VX** 
20ELTA*=  A«F7.6*1X*A2*/*16X*ATHETA>  A ,F7 .6* 1X*A 2) 

107  FORMATtSFlO.O) 

lOB  FORMAT (/*BX*AEXPERIMENTAL  VELOCITY  PROFIUA*/,  19X*  I2*A  POINTSA*/*! 

1<»X*AV  (A,A2*A)  A*6X*AU/UEA) 

109  F0RMAT(13X*F7.6,5X*F6.6> 

110  F0RMAT»//,5X,AAFTER  A,I2,A  ITERATIONS  IN  ORSRCH,  FlX)a  A*E13.6«A  k 
CITM  AA*I2,2X*APARAHETER  FITA) 

111  FORHAT(5X*ACAPK  «A *E13.6*2X* AS  «A* E13. 6* 2X* AKAPPA  «A*E13.6*2X* AUTA 
CU  sA*El3.6*/l 

112  FORMAT (5X*ACIN«  A*Et3 .6«5X *AT0PLUS«  **E13.6I 
113  F0RMAT(/*BX«A8ETAC  «  A*El3*6) 

116  FORMAT  (//*SX*A(1ATCH  CONBITION  CRITERIA  FOR  VELKITY  PROFILEA) 

115  F0RMAT{5X*AETAM  »  A*F7.6«6X*AN0  »  A  * l6*6Xt AETAO  «  A*F7.6,) 

116  FORMAT  (1H1*13X*  66M*  EXPERIMENTAL  PROFILE  •  ANALYTICAL  PROFILE 

1  •  U>  DEVIATION  •*/*16X«  1N**26**  1H**22X,  tH«*16X*  lN*f/*( 

2X*  2HY**6X*  1H»*5X*  5HUD/UE*7X*  3HU06*6X*  lH*t9X*  6HUA/UE*5» 

3*  3HUA»*6X«  1N**3X*  9HU0^  >  UA»*6X*  IH*) 

117  FORMAT  (5X*F7.2*2X*  IH* *SX*F7.6*6X*F7.2*3X*  IN** 3X*F7.6*2X«Fr.2* 

13X*  1H**3X*F9.6*6X*  1H«) 

118  FORNATI1H1*////////*23X«5A10*/,61X*A(A*A7*A  UNITSIA* 

•  //*//*39X*A£XPERINENTAL  VALUESA*7*  ltX*AIDA*  6X*AXSTAA*5X*AUEA«5)i 


»,«UTAU«t8X««NU*,6X.«eCT««,6X,«0ELr(«,4X,«0ELTA •A«3Xt*THETAX,SX, 
♦ATUAI 

119  FORHATI10XtI4t3X«F6.3»2X»F6.2t3X.F6.4«3X.FB.7t3X.F6.3t4l2X*F7.4n 

120  FORHAT  (////« 15X«#tJNSTEA0Y  MALL  ♦  SIMILARITY  HOCEL  -  FULL  FROF 

♦ILEF.I4.2X.FPARAHETER  FZTF,//«llXt  2HIO*5X«  4HBETA«SX«  4HUTAU« 
46Xt  lHS«aX.  IHK.aX*  SHKAPPA,2X»  7HEPSlLON«liX.  3HT0»t8X» 

42HCII 

121  FORMAT  (lOX«I4«lX«F8.4t2X,F6.3.2X*F7.3.2XtF8.6t3X«F7.4,2(lX,F8.61 t 
13X,F6.31 

122  F0RMAT(/«2ax,FMEAN  EPSILON  OVERF»I 4» 2X . FSTATIONS  «««F8.6I 
130  FORNATIFIO.O) 

200  F0RMAr(//«9K«FLEN6TH  SCALESFl 

201  FORMAT  (SXtFOELIN  s  * ,F7. 5,a>X«F0EL3UT  »  F,F7.SI 
RETURN 


oooooooooooooooooo 


SUBROUTINE  ORSRCH  f N, NOIN« F, F4,FS« X,  OX, EPSI ,MA X,EPSIF,ITER t IPRINT* 
HER) 

EXTERNAL  F 

OINENSION  FSCNOIN),  x:NOIN),  OXtNOIN),  EPSICNOIN) 

ORSRCH  -  DIRECT  SEARCH  ROUTINE 
VARIABLES 

N  -  NUHBER  OF  PARAMETERS 
NOIN  •  OIHENSION  OF  ARRAYS 
F  -  FUNCTION  NAME 
Fl»  -  FINAL  VALUE  OF  F 
FS  -  VECTOR  OF  INTERMEDIATE  F  VALUES 
X  -  VECTOR  OF  INITIAL  AND  FIVAL  PARAMETER  VALUES 
OX  -  VECTOR  OF  INITIAL  STEP  SUES 
EPSI  >  VECTOR  OF  CONVERGENCE  :?ITEREON 
MAX  -  MAXIMUM  NUMBER  OF  ITERATIONS 
EPSIF  -  FINAL  CONVERGENCE  CRITETEAN 
ITER  •  NUMBER  OF  m RATIONS 
IPRINT  -  PRINT  CONTROL 
lER  •  ERROR  FLAG 

IER>0 

ITER«IER 

IF  (IPRINT.EQ.l)  MRITE  <6,120) 

Fi»«F(X) 

101  IsQ 

102  I>I»1 
MSTEPaO 
X3»X{I) 

F3*F4 

IF  (IPRINT.EQ.l)  MRITE  <6,121)  IT- R,FT, (X( J) ,0X(J) ,Jal,NI 

103  ITER=ITER»1 

IF  (ITER. GT. MAX)  GO  TO  110 

NSTEP»NSTEP»1 

X2aX3 

F2aF3 

XtI)aX(I)«OX(I) 

X3aX(I) 

F3aF<XI 

IF  (IPRINT.EQ.l)  MRITE  <6,121)  ITE?,F3, (X ( J) ,OX(J) ,Jal,N) 

IF  (F3>F2)  104,104,106 

104  XlsX2 
F1*F2 

GO  TO  103 

106  IF  (NSTEP-1)  106,106,107 

106  XlaX3 
F1»F3 

OX<I)«-OX(I) 

X{I)«X{I)»OX(I) 

X3«X2 
F3«F2 
GO  TO  103 

107  ITER«ITER*1 

IF  (ITER.GT.MAX)  GO  TO  110 

X22«X2*X2 

X32«X3»X3 
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F2HF3»F2-F3 

AFAC=F1*«X2-X3I-X1*F2MF3»F2*X3-F3**2 
IF  {AFAC.EQ.O.t  60  TO  119 

8FAC*X1*X1*F2MF3-F1»IX22-X32>>X22»F3-X32»F2 

X4a-BFAC/(2.*AFACI 

XdlsXb 

FI»aF{Xl 

IF  (IPKINT.EQ. 1)  HOITE  16.122)  ITER. F6« fX (J ) «OX(J> «J>1,N> 
0FaFii-F2 

IF  <ABS(2.«0F/(F4>F2)).LT.EPSIF>  SO  TO  114 
0ELXaX4‘‘X2 

IF  IA6S(0ELX).LT.EPSI<I>)  GO  TO  114 
IF  (OXdl.LT.O.)  OCLXa-OELX 
IF  (OF)  106.10a.lll 

108  IF  (OELX)  109.109.110 

109  X3sX2 
F3=F2 
X2aX4 
F2aF4 

60  TO  107 

110  XlsX2 
FlaF2 
X2*X4 
F2aF4 

60  TO  107 

111  IF  (OELX)  112.112.113 

112  XlaX4 
FlaF4 

GO  TO  107 

113  X3aX4 
F3aF4 

GO  TO  107 

114  IF  (H.EO.l)  RETURN 
FSd)aF4 

IF  (I'D  102.102.119 

115  IF  d>N)  102.116.116 

116  IF  (ABS(2.*(FSd)-FSd-l)l/(FSd)»FSd-l)))  .LT.EPSIF)  RETURN 
00  117  Jal.N 

117  ox:j)aoxu)/i.i 
60  TO  101 

118  lERal 
MRITE(6.124)  NAX 
RETURN 

119  IERa2 

NRITE  (6.123)  XI ,Fl.X2*F2.X3.F3 
RETURN 

120  FORNAT  (1H1.3X.  5HCTCLE.8X.  1HF.13X.  4HXd).9X.  9H0X(I)) 

121  FORNAT  I  9N  E  .I3»3X.E13.6.6(2X.E13.6)) 

122  FORNAT  (  9H  I  « I3.3X.E13.6.6(2X. E13.6) ) 

123  FORNAT  (9Xt  4HX1«  .E13.6*  6N  Fla  .E13.6.  6H  X2a  «E13.6. 
1F2«  .E13.6.  6N  X3a  .£13.6.  6H  F3a  .£13. 6) 

124  FORNAT (5X.AORSRCH  FAILED  TO  CONVERSE  IN  A. I4 .8ITERATI0NS8) 

END 
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FUNCTION  F:X» 

DIMENSION  X(3) 

C0HN0N/VAR/IVAR(3>  *X0 131 

COHNON/OATA/NPAR«NOF«XNUtUI*OUIOX«r3(90l«UO(90ltOLSTR*UA(90) » 
»YOP(.:ioQoi  «uAPt.(t«oai  «nppt 

CONMON/UOUT/NO .ETA  C90  01 . OETA  UOO) * %00  >  *ETAN. ETAO .COUT.CIN. 
IS.TNOTP.OELIN.OELOUT.USrAR.BETAC 
DATA  EPS/1 .£-10/ 

C  FUNCTION  F  COMPUTES  THE  ROOT-NEAN-SaUARE  ERROR 
C  BETHEEN  THE  MEASURED  EXPERIMENTAL  VEL3CITT  PROFILE  DATA 
C  AND  THE  THEORETICAL  VELOCITY  PROFILE.  THE  THEORETICAL 

C  PROFILE  CONSISTS  OF  THE  UNSTEADY  MALL  LAYER  MODEL 

C  FOR  THE  INNER  LAYER  AND  A  SELF-SIMILAR  PROFILE  MODEL 
C  FOR  THE  OUTER  REGION. 

C 

C  ARRANGE  VARIABLES  INCOMING  FROM  DSRCM  IN  CORRECT  ORDER 
C 

JCsl 

DO  10  1*1.3 

IFdVARlII.EQ.OlGO  TO  It 
X0{I)*X ( JCI 
JCsJC»l 
10  CONTINUE 

CALL  UTAUF  TO  DETERMINE  UTAU 

XKAP«X0(31 
S*X0(2I 
CAPKxXOdI 

CALL  UTAUF (CAPK.XKAP) 

COMPUTE  RMS  ERROR 

JLOC«2 
F=0. 

UASlIxO.OOO 
DO  20  1*2. NOP 
VOP*TO(II/OELIN 
ETAO*VO(I>/OELOUr 
C 

C  COMPUTE  THE  DEFECT  PROFILE  AT  THE  DATA  POINT  EIT»€R 
C  FROM  THE  SERIES  SOLUTION  FOR  ETAD  LESS  THAN  ETAN 
C  OR  BY  INTERPOLATION  OF  THE  NUMERICAL  SOLUTION  FOR 
C  ETA  GREATER  THAN  ETAN 
C 

IFtETAO.GT.ETAmGO  TO  3t 

CALL  UlSER (XKAP.BETAC.COUT.ETAO.U1.U10.EPS) 

HAKE«U1 
CO  TO  AO 

30  DO  SO  JxJLOC.NO 
Xl«ETA(Ul-ETAO 
X2*ETA(J-ll-ETAO 
X3«X1*X2 

IFIX3.lt. O.IGO  TO  60 
SO  CONTINUE 
MAKE*0. 

GO  TO  AO 
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60  JLOC-J 

HAKE>U:jl-tu:jl-U(J*lll*Xl/OET«(J>ll 
60  MAKEaMAKE-ALOGtETAOI/XKAP-COUT 

UAiI>aUSTAR«<UP(VOP*StTNOTP.ClN.X<AP,0. I»NAKE) 
rT3UO(II-USTAR«(UP(VOP*S*TNOTP»CIN*XKAP,0.)  «’MAIO 
20  F»F*FT*FT 

F*F/FtOAT«NOPI 
70  CONTINUE 
RETURN 
END 
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FUNCTION  EOVVIS(XKAP*CAFK«ETAI 
DATA  N.XN/<»,a.25/ 

FUNCTION  EOYVIS  EVALUATES  THE  EDOV  VISCOSITY  FUNCTION 
OR  ITS  DERIVATIVE  FOR  THE  OUTER  LATER  OF  A  TURBUIENT 
BOUNDARY  LATER.  THE  FUNCTION  APPROACHES  KAPPA»ETA 
FOR  SHALL  ETA,  EXPONENTIALLY  QUICKLY.  FOR  LARGE  ETA 
THE  EOOY  VISCOSITY  FUNCTION  APPROACHES  THE  OUTER  CONSTANT 
K  ALGEBRAICLY 

INPUT  PARAHETERS  XKAP  •  VON  KARHAN  ACONSTANTF 

CAPK  •  OUTER  REGION  ACONSTANTA  K 
ETA  >  SCALED  OUTER  REGION  COORDINATE 

I^O 

GO  TO  10 
ENTRY  EDYVISO 
I«1 

10  C>ICAPK/XKAP  )**N 
Xl*EXPt-C/ETA**NI 
X23l.>Xl 

IF(I.EQ.O>EOYVIS«CAPK 
IFd.EQ.llEOYVISaO.OO 
IFCETA.6T.10.0I  GO  TO  20 
IF ( I .EQ. 01 EOYV IS»XKAP*ET A*X2**  XN 
IF(I.eQ.llE0YVIS«XKAP*X2**XN*U.-;*Xl/<X2*ETA**Nn 
20  RETURN 
END 
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SUBROUTINE  UOUTERi XKAP,CAPK1 

OINENSION  UC(<ta0«2l«COIZI«8A(2),F(l>00)*DELI<»00  l«  A  (%00)  «8 1 4C0I  t 
1014001 

CONHON/UOUT/NO*ETA(4  00I,OETA(400I  .Ut400)«ETAN,ETA0«COUTtCINt 
lS,TNOTPtOELlN«OELOUT«USrAR,BETAC 
DATA  EPS0<ALPHA«EPS/l.E>10«l.a2.1.£*ia<' 

SUBROUTINE  UOUTER  COMPUTES  A  NUHERIIAL  SOLUTION  FOR  THE 
OUTER  REGION  VELOCITY  DEFECT  PROFILE  GIVEN  THE  IfTUT 
VARIABLES  XKAP  -  VON  KARHAN  FCONSTANTF 

CAPK  >  OUTER  REGION  EOOV  VISCOSITY  ACONSTANTF 
BETAC  -  CLAUSER  PRESSURE  GRADIENT  PARAMETER 
NO  •  NUMBER  OF  MESH  POINTS  8ETHEEN  THE  HATCH 

MITH  THE  SERIES  SOLUTION  (FOR  SMALL  ETA  I  AND 
THE  BOUNDARY  LAYER  EDGE. 

ON  INPUT 

ETA (I)  -  MESH  POINTS  BETWEEN  ETAH  AND  ETAO 
DETAdl  -  variable  STEP  SIZE  (INITALLY  UNIFORM  AND 
EQUAL  TO  H  FOR  FIRST  6  NESH  POINTS*  THEN 
SUCCESSIVELY  INCREASING  BY  A  FACTOR  ALPHA) 

U(II  >  NUMERICAL  SOLUTION  FOR  VELOCITY  DEFECT  AT  ETA (It 
ETAM  •  MATCH  POINT  BETWEEN  SERIES  AND  NUMERICAL  SOLUTION 
ETAO  •  LARGEST  VALUE  OF  EIAd) 

COUT  -  OUTER  REGION  LOG'LAV  FCONSTANTF 

FIX  THE  MATCH  POINT 

ETAM«0.5*CAPK/XKAP 

FIX  THE  OUTER  VALUE  OF  ETA 

X1>1.»2. •BETAC 
X2-X1-1. 

X33-(2.*BETAC»1.I/(2.«BETAC) 

IF(X1.GT.O.IETA03SQRT(-2.«CAPK*AL3G(EPSO)/X1) 

IF (Xl.EQ. 0 .1 ET AOa-SQRT (ABS (CAPK/XZ)  I  • ALOGIEPSO ) 
IF(X1.LT.O.IETAO>SQRT(ABSICAPK/X1>I»EPSO«»(X3) 

define  the  mesh 

NlaNO-1 

ETAdlsETAN 

H»tETA0-ETAN)F :4.»(1.-ALPHA**(N0-5)  ) F (1 .-ALPHA!) 

DO  10  1*1. N1 

IF(I.LT.6)0ETA(II«H 

IF(I.GE.6IOETA(I)«ALPHA*DErA(I-l} 

eta;i»i?seta:i)»oeta(I) 

10  CONTINUE 

CALCULATE  THE  ELEMENTS  OF  THE  TPIDIAGONAL  MATRIX  FOR 
THE  NUMERICAL  SOLUTION  FOR  THE  DEFEST  PROFILE 

XP*ETAIH»0.5*H 

EPaEOYVIS(XKAP,CAPK*XP) 

EPOsEDYVISDIXKAP.CAPK.XP) 

P JP« ( EPD •X 1* XP ! / EP 
QJP*Y2/EP 
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00  20  Ja2«Nl 
GAH»t. 

IFiJ.GE.eiGAHsALPHA 

X3s0.2G*OETA:j1«OETA{J1 

Xi»=Q.5*0ETA(JI 

XMaXP 

XP=£rA<JI»X4 

EPsEOVVIS(XKAPtCAPK*XPl 

EP0-E0VVIS0(XKAP*CAPK«XP) 

PJM*PJP 

PJP=(EP0*X1»XPI/EP 

QJHsQJP/GAM 

QJP*X2/EP 

A(J|3-1.-6AN-X4*(PJP>PJNI»X3>IQJP»9JN) 

8(JI=1.*X4*PJP ♦XS^OJP 
20  C(J>=GAH-X4*PJM*X3*QJM 

DEFINE  TMO  INITIAL  GUESSES  FOP  GOUT 

C0(1I«1. 

C0t2lsa.75»C0(l) 

CALCULATE  TMO  NUMERICAL  SOLUTIONS  FOR  U  FOR  ETA  GREATER 
THAN  ETAN  FOR  GUESSED  VALUES  OF  COUTER 

00  30  K3l«2 

CA LL  U1 SER ( XKAPt  BETAC  tCO (Kl , ET AH*  JC 1 1 1 K ) «  UCOS*  EPS1 
F(1I>0. 

0EL(l>«UCa*K) 

00  40  J32«N1 

xi*a:ji*c:ji*fij-i» 

F(J>=-B(J>/X1 

40  0EL(JI«-CIJI*0EL(J>1)/X1 

BACK  SU8STIT10N 

UC(N0«K)«0. 

00  90  Ja2«Nt 
J1=N0-J»1 

uc:ji«k|3F(jii«uc(ji»i*k)«’OELiji> 

50  CONTINUE 

CALCULATE  DERIVATIVE  OF  NUMERICAL  SOLUTION  AT 
CTA>ETAN  HITH  SLOPING  DIFFERENCE  FORMULA 

UCON-(-274.«UC(l>KI*600.*UCI2*K>>900.*UC(3«K) 
1»400.*UC(4*KI-1S0.*UCI5«K) *24.*UCI6*K)  I  /  :120.»HI 
CALCULATE  DIFFERENCE  DETICEN  NUMERICAL  AND  SERIES  DERIVATIVE 
SOLUTION  AT  ETA«ETAN. 

8A(KtaUC0S>UC0N 
30  CONTINUE 

C0H8INE  TMO  NUMERICAL  SOLUTIONS  TO  OSTAIN  THE  TRUE 
VALUES  OF  Util  AND  COUT 

8l«-8A{2l/:BAtll*0Af2ll 


82=1.-81 

C0UT=ai«C0(l)»B2«C0(2t 
00  60  J=1«N0 

U(Ji*ai«UC(Jtll*82*UC(J«2) 
60  COHTINUE 
RETURN 
END 
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SUBROUTINE  UiSER(XKA(>«BETAC«COUT,£Tt,Ul,UtOtEPSI 
DATA  NHAX/IOO/ 

subroutine  UlSER  EVALUATES  THE  WELOSITY  DEFECT  FUNCTION  Ul 
ANO  ITS  DERIVATIVE  UtO  AT  THE  INPUT  VALUE  OF  ETA,  FRON 
A  SERIES  SOLUTION  lOBTAlNED  BY  THE  NETHOO  OF  FROBENIOUS) 

WHICH  IS  VALID  FOR  SHALL  ETA  ANO  HHISH  SATISFIES  THE 
ordinary  DIFFERENTIAL  EQUATION  FOR  THE  SIMILARITY  DEFECT 
PROFILE  Ul. 

INPUT  PARANETERSl  XKAP  >  VON  KARNAH  ACONSTANTF 

BETAC  >  CLAUSER  PRESSURE  GRADIENT  PARAHETER 
COUT  -  OUTER  REGION  LOG-LAM  « CONSTANTA 

EPS  -  convergence  tolerance 

UlsUlO^O. 

Xt3ALOG(ETAI 

X2a2.*BETAC 

X3»X2*1. 

FACAal./XKAP 

FACa«COUT 

ASUH>FACA 

BSUN>FAC8 

ASUNOaBSUNOaO. 

DO  10  IsltNNAX 
XA-FLOAT(I-l> 

X5sXA»l. 

X6*1./JX5*X5*XKAPI 

xra(X2»XA«X3l*X6 

FACB«(-FACB«X7»<X7-X6I«FACA/X5I«ETA 
FACA«-FACA*eTA*X7 
ASUNaASUN«FACA 
BSUN«8SUH»FAC8 
TESTAaASUND 
TESTBaSSUNO 
ASUNOaASUNO»X9*FACA 
SSUHOaeSUHO»XS«FAC8 
IF(I.LE.2)GO  TO  10 
IFtASUND.EQ.O.QlGO  TO  5 
IF(A8SI1.-TESTA/ASUN0|,GT.EPSIG0  to  10 
9  IF(ABS(1.-TESTB/BSUN0I.LT.EPS)G0  to  20 
10  CONTINUE 

PRINT  30,  NNAX, XKAP, BETAC, COUT, ETA, EPS 
30  FORNATI//,lX,«SERIES  CALCULATION  IN  UlSER  HAS  NOT  CONVERGED  AFTERP 
»I<i,2X,«TERNS«t//tlX,PXKAPa  «,E12.G,3X ,ABETACa  P.E12.S, 

*3X,,«COUTa  A,El2.5,3X,«T0LERANCEa  A,E12.9,//) 

STOP 

20  UlaBSUN«Xl*ASUH 

U10a{BSUNO*Kl*ASUNO»ASUN)/ETA 

RETURN 

END 
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SUBROUTINE  UTAUF(C«PK*XKAP) 

COHHON/OATA/NPAR«NOP«XNU*UI.OUIOX,r3(90)*UOI90l«OLSTR«UA(90) , 
♦YOPLClOaO) tUAPLtlOOOl «NPPT 

COMHON/UOUT/NO«ETA (VOO)«OETA(4aOt . Jt400)»ETAH,ETA0«COUT«CIN, 
IStTNOTP.OELZN.OELOUT.USTAR.BETAC 
DATA  EPS.NAX/1^-8,20/ 

UTAUF  COMPUTES  THE  VALUE  OF  UTAU  ITERATIVELY  USING 
THE  MATCH  CONDITION  AND  ALSO  THE  OUTER  REGION  DEFECT 
PROFILE 

EPS  -  EXIT  TOLERANCE  FOR  TMO  SUCCESSIVE  ITERARES  FOR  UTAU 
HAX  -  MAXIMUM  NUMBER  OF  ITERATIONS  IV  UTAUF 

CINsCINF(S«TNOTP,XKAP.O.*l.) 

REXsUI/XNU 

USTRA^USTAR 

OO  10  Isl,MAX 

OE  LOUTa  OLSTR/USTRA 

0ELIN>1./(REX«USTR At 

BETAO-OUIOX«OELOUT/IUI«USTRA> 

IFIBETAC.LT.'.SOI  BETACs'.SO 

CALL  UOUTERtXKAP.CAPKI 

XAal./ ( ALOG(OELOUT/OELINI /XKAP»CIN-COUT) 

IFd.NE.llGO  TO  20 
FB*USTRA-XA 
USTRB«USTRA 
USTRAaXA 
GO  TO  10 
20  FAaUSTRA'XA 

XlaUSTRA-FA*(OSTRB-USTR»l/IFB-FA» 

TESTsABSd.-USTRA/Xll 
IFITEST.LT.EPSIGO  TO  30 
IFtABSfFBI .LT.ABS(FAI)  GO  TO  15 
USTRB-USTRA 
FB=FA 

15  USTRA«X1 
10  CONTINUE 

MRITE(6t40)  MAX«USTRA,X1 

40  FORMAT(/«lX«XITERATION  HAS  FAILED  T3  CONVERGE  IN  UTAUF« 

1AFTERP«I4«»  ITERATIONSP*/«lX»ALAST  TWO  ITERATES  FOR  UTAU  AREF« 
22E16.af/) 

30  USTARsUSTRA 
RETURN 
END 
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C 

C  *  •  •  FUNCTION  CINF  USES  THE  NALL  CONPATIBILITY  CONDITIONS  TO 
C  *  •  •  CALCULATE  THE  INNER  REGION  CONSTANTS  CINER  AND  TNOTP  FOR  A 
C  •  •  •  GIVEN  VALUE  OP  S  FOR  BOTH  THE  TEHPERATURE  AND  THE  VELOCITY 
C  •  «  «  PROFILE  IN  THE  INNER  LAYER  OF  A  TURBULENT  BOUNDARY  LAYER. 

C  •  •  •  ROOTPR  a  i.  FOR  THE  VELOCITY  PROFILE  CASE. 

C  •  •  •  PPLUS  =  0.  FOR  THE  TEMPERATURE  PROFILE  CASE. 

C 

DATA  GM,ROOTP1«EPSI«ITMX/>0. HO 45393 V8109179*!. 772  46386690 952. l.E-1 

*0.20/ 

JCaO 

Xl=0.5*S*ROOTPR»ROOTPI*KKAP 

S2*S»S 

X2a2.*XKAP*S2*PPLUS/3. 

ALa£XPJ-l.-Xl-X2> 

00  10  Kal.ITHX 

XSaALOGCALI 

X4aAL»l. 

AL2aAL*AL 

X5aSQRT(l.-AL2) 

Ca  AL 

ALaAL-<Xl*X5«^X4*X3n.-AL  ♦X2*IAL2*4.*ALH.»  7X4» /I-AL*X1/X5» 

*  X3«l./AL»X2*tAL2»2.»AL»3.)/(X4*X4> ) 

IF(AeS(AL/C>l.l*LE.EPSXI60  TO  20 

10  CONTINUE 
JOl 

20  AL2aAL«AL 
X4al.-AL2 
TN0TPaAL2*S2/X4 
X6aSQRT(X4l 

CINF«X1/X5 *6N-0,5*ALOGIS2»TNOTPI  »1.F(AL«’1.I 

♦  ♦I0.25*X2*«1,-9.*AL2» -PPLUS*TNOTP*XK APa :3.*  «.2-4.»AL»l .> /3. 1/X4 
CINF»CINF/XKAP 

IFCJC.EQ.OIGO  TO  30 
HRITE(6,40IK.TNOTP,CINF» AL.C 
STOP 
C 

40  FORNATdX.  3QHNO  CONVERGENCE  IN  ;iVF  AFTER.IH.  IIH  ITERATIONS./. 
»1X,  7HTN0TP  a.E15.5.10X.6HCINF  «.E15. 5. lOX . 4HAL  «.E15.5.1QX.5HALP 
4a.E15.5./J 
C 

30  RETURN 
END 
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c 

c  • 
c  • 
c  • 
c  • 
c  • 
c  • 
c  • 
c  • 
c  • 
c  • 
c  • 
c  • 
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FUNCTION  UP(VF,S«TN0TP*CIN*XKAP.P3VUS1 

•  *  FUNCTION  UP  CALCULATES  THE  TINE-HEAN  TEMPERATURE  OR  VELOCITY 

•  •  PROFILE  IN  THE  MALL  LAYER  OF  A  TURBULENT  BOUNDARY  LAYER 

•  •  FOR  SPECIFIED  VALUES  OF  THE  PARAMETERS  LISTED  BELOM. 

•  •  YP  -  Y«« SCALED  HALL  LAYER  COORDINATE 

»  •  S  -  S  .CYCLE  TlHE  PARAMETER 

•  •  TNOTP  -  T0*,SINILARITV  PARAMETER 

•  »  CIN  -  INNER  REGION  LOG'LAW  CONSTANT  Cl  OR  Bl 

•  •  XKAP  -  KAPPA  (VON-KARMAN  tCONSTANTAI  OR  KAPPA-THETA 

•  •  PPLUS  -  HALL  LAVER  PRESSURE  GRADIENT  PARAMETER.  (PPLUS 

•  •  IS  lERO  FOR  THE  TEMPERATURE  PROFILE  CASE.» 

•  •NOTElA  CALL  TO  UP  SHOULD  NORMALLY  BE  PRECEDED  BY  A  CALL  TO 

•  •  CINF  HHICH  COMPUTES  CIN  AND  TNOTP.  GIVEN  S. 

DATA  X1.SRPI/-0.40A53934B1091T9.1.TT2A53BB0905S16/ 

P  t  XI a-2 . ’PPLUS* I X*TNOTP» /3 , 

R ( XI aCIN* ( 0. 9*  ALOG (X  *TN0TPI-X1 1 /X<  AP *0 . 5*PPLUS*tS2  *2 .•TNOTPI 
(i<X.Y,ZAI»<2.*X*X*l.l*r*2.*X*2A/SRPI 

I(X.Y,ZAI»A.*< (2.*X*X*l.l*Xi:x)*X*XIP:xl-0.125*SRPI*tG.*X*X*l.l*Y 
♦  -0.75*X*ZAI/ISHPI*XKAP» 

HtX.Y,2AI»IX**%*3.*X*X*0,75)*r»X*IX*X’2.5)*2A/SRPI -3.*X*X 


•  •  *  PRECIS  IS  THE  VALUE  OF  X  SUCH  THAT  EXP:-X*XI  MAV  BE  COMPUTED 

•  •  *  MITHOUT  INCURRING  AM  UNBEFLOH. 


PRECI$*29.936tA99 

S2*S*S 

TPS2»S2’TN0TP 
M»0.5*YP/SQRT<S2 ’TNOTPI 
H0s(l.5*rP/SQRT  JTNOTPI 
ERFH«CRF{H1 
ERFH0«ERF(HQI 
EXPH«EXPH030. 

IF<H,LT.PR£CISI£KPHaEXPf-H*MI 
IFiHO.LT.PRECISI  EXPHO*EXPA-HO*HO> 

UP=TPS2*  <R <S2I *0 (H.ERFM.EXPHI *2 IH, • RFH.EXPMI I 
’  -TNOTP* (R ( 0 .1  *0 (HO .ERFHO.EXPHOI  *Z (HO .ERFHO. EXPHOI I 
IFCPPLUS.EO.Q.IGO  TO  10 
IFCH.LT.PRECISIGO  TO  20 
UP«UP-fl . 5*  CS2*S2  *2 .*S2*TN0TPI*PPLUS 
GO  TO  LO 

20  IFIHO.LT.PRECtStGO  TO  30 
MH»H (H.ERFM.EXPHI 
HHfl»H0**A*.75 
60  TO  bO 

30  HHsHtH.ERFH.EXPHI 

HH0«N<H0. ERFHO.EXPHOI 
A 0  UP*UP*TPS2*P (S2I •HH-TNOTP*P 1 0. 1 •MHO 
10  UP«UP/S2 
RETURN 
ENO 
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FUNCTION  XIIX) 

DATA  SRPItGANa«CPSI/i.772%53850905516«0.57721566490153286«l.E-iO/ 
C 

C  •  •  •  FUNCTION  XI  EVALUATES  THE  TRIPLE  INTEGRAL  IN  THE  UNSTEADY 
C  *  •  •  HALL  LAVER  MODEL. THE  TOLERANCE  EPS I  IS  THE  NUHBER  OF 
C  •  •  •  SIGNIFICANT  FIGURES  DESIRED  FOR  XI  AND  ITS  DERIVATIVE 
C  *  •  *  XIP  AND  IS  MACHINE  DEPENDENT. 

C 

XI=0SIF»1 

IFtX.LE.O.)  RETURN 
GO  TO  1 
ENTRY  XIP 
XI«0.SIF32 
IFtX.LE.O. I  RETURN 
1  X2»X*X 

FAC»2.»X2 

M>taa 

SUM«fl. 

SUHTsO. 

TERN«1. 

IF(X.GE.5.30I  GO  TO  110 

if:if.eq.2)  go  to  lOO 

TERMS X 
ALPHAsl. 

DO  2  Isl*H 

TERM»TERH*FAC/FL0AT(2*1P1I 
S  UMs  SUM  ♦  TE  RH*  A  LP  HA 

IF(AeS((SUN-SUNT)/SUN).Lr.EPSI>  60  TO  3 
ALPHAsALPHA»1. /FLOAT  I !♦! ) 

SUMTsSUN 

XIa0.2S*EXP<-X2J*SUM 

RETURN 

00  DO  h  Isl.H 

TERMsTERM* FAC/FL OAT t 2»l- 11 
SUM>SUM»TERH/FLOAT  CII 

IF(ABS((SUN-SUMTI/SUHI.LT.EPSI>  GO  TO  S 
SUMTaSUM 

X1*0.25*EXP(-X2I*SUN 

RETURN 

110  IF(IF.EQ.2I  GO  TO  120 

DO  6  laltM 

TERMS  TERM*  FL0AT12*I-H /FAC 
TERMAsTERM/FLOATCI) 

IF(TERNA.LT.EPSI)  GO  TO  7 

6  SUMsSUM^TERMA 

7  XIsSRPI«(AL0G(X2)^GAH0>SUN)/8. 

RETURN 

120  00  8  Isl,H 

TERNaTERM*FLO AT  C  2» I- 1» /F  AC 
IFITERM.LT.EPSII  60  TO  9 

8  SUMsSUM^TERM 

9  XIsSRPI*(l.»SUMl/(l».»X> 

RETURN 

END 
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SUBROUTINC  LOGAXINCYC.XLC.HGHT) 

C  SUBROUTINE  LOGAX  PLOTS  THE  X>AXIS  L3GARITHHICALL V 
C 

C  INPUT  PARAMETERS 
C  NCYC  s  NUMBER  OF  CYCLES 
C  XLC  3  X  LENGTH  OF  ONE  CYCLE 

C  H6HT  «  HEIGHT  OF  AXIS  LABEL 

C  30  FOR  NO  LABEL 

DIMENSION  XLOGiAl 
TCKlsHGHT 

IFtHGHT.LE.a.t  TCXt«XLC/EO. 

TCIC2aTCICl/2. 

C  SET  UP  EIGHT  TIC  HARKS  PER  CYCLE 
DO  1  l3t,S 

1  XLOG;i)«ALOG10tFLOAT(I»l>)*XLC 
DO  3  NxltNCYC 

XNsN 

XNHlsXN-l. 

XLCNH1«XLC*XNM1 

IF(N.EQ.l)  CALL  PLOT (XLCNMl.TCKl.S) 

CALL  PL0T(XLCNHl«0.*2t 
00  2  l3l«a 
X«XLCNM1»XL0G(I} 

CALL  PLOT<X*0.,2) 

CALL  PL0T(XtTCK2*2) 

2  CALL  PLOT<X«0.«2I 
XNLC«XN*XLC 

CALL  PLOT(XNLC«0.«2) 

3  CALL  PL0TIXNLC«TCX1«2» 

IFtHGHT.LE.d.l  GO  TO  5 

NUM»10**NCYC 

FACs3.*HGHT/7. 

TNUN« FLOAT IS-NCYCl  ’FAC 

NCYCP1=NCYC*1 

00  A  N«t«NCYCPl 

NNUM-NCYCPl'N 

XNUNaFLOAT (NMUMI 'XLC-TNUN 

TNUNaTNUN^FAC 

YNUM«-1.5»HGHT 

CALL  NUN8ER(XNUN«YNUH»H6HT«NUHtO.»2HI5) 

A  NUNbNUM/10 

XSYM»XLC»FLO AT (NCYCI /2.-2.»FAC 
YSYM»-3.*HGHT 

CALL  SYMBOL ( XS YN*rSTN, HGHT «2HY »f 8. 1 2 > 

5  CALL  PLOT(0.fO.«3l 

RETURN 
ENO 
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SUBROUTINE  SHFTVAX  (NLtNRl.TPPtXL*I>«CPL«HGHT) 

C  SUBROUTINE  SHFTYAX  PLOTS  THE  Y-AXIS  FOR  A  MULTIPLE  NUMBER 
C  OF  CURVES  ON  THE  SAME  PLOT  MITH  SHIFTED  Y-AXIS  ORIGINS 
C 

C  INPUT  PARAMETERS 

C  NL  >  NUMBER  OF  POINTS  LABELED  NON-IERO 

C  NPLTPP  >  NUMBER  OF  POINTS  LABELED  ZERO 
C  (EQUAL  TO  NUMBER  OF  AXIS  SHIFTSl 

C  XL  3  LENGTH  BETNEEN  AXIS  TICS  (SPACING) 

C  INCPL  >  a 

C  HGHT  s  HEIGHT  OF  AXIS  LABEL 

TCKaHGHT 

IFfHGHT.LE.O.)  TCKaXL/m. 

N«NPLTPP^HL-1 
CALL  PL0T(TCK«a.«3) 

CALL  PLOT(O..O.tn 
YaO. 

DO  1  Ial«N 
YaY»XL 

CALL  PL0T(a.tY«2) 

CALL  PL0T(TCK«Y«2) 

1  CALL  PLaTtO.*Y.2l 
IFfHGHT.LE.O.)  60  TO  % 

NUNaNL* INCPL 
FACaG.*H6HT/7. 

XNUN«-2.5«FAC 

YNUNay-H6HT/2. 

00  2  I«1,NL 

IF(NUN.LE.9B)  CALL  NUNBER(XNUH*YNt)Nf  HGHT»NUN»0  .•2HI2) 
IFtNUM.GT.BB)  CALL  NUMBER(XNUH-FA3«TNUN»HGHT.NUM«0.*2HI3) 
YNUNsYNUM-XL 

2  NUHa NUN- INCPL 
XNUM=-l,5*FAC 
00  3  I«l, NPLTPP 

CALL  NUN8ER(XNUN*YNUN«H6HT«NUH«0.t2HIl) 

3  YNUMaYNUM-XL 
XSYM=-3.*HGHT 
YSYM*XL*FL0AT(HI/2.-FAC 

CALL  SYMBOL(XSTN»VSYN«HGHTt2HU»«90.t2> 

I*  CALL  PLOT(a.tO.»3) 

RETURN 

ENO 
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APPENDIX  D 
TEST  CASE 

The  turbulent  boundary  layer  prediction  code  presented  in 
Appendix  C  is  set  up  to  guide  the  user  in  understanding  the 
code's  operation.  Connents  are  used  to  explain  input  variables 
and  to  denote  where  major  operations  are  taking  place. 

The  user  supplied  code  input  begins  with  "PROGRAM  MAIN"; 
an  integer  value  for  each  of  the  model  parameters  K,  S,  and  k 
is  read  in  a  312  format.  This  integer  value  determines  if  the 
parameter  is  to  be  optimized  or  not;  an  integer  value  of  0  indi¬ 
cates  no  optimization  for  the  corresponding  variable  and  1  indi¬ 
cates  parameter  optimization.  The  next  variables  to  be  read 
are  "ITERM",  "IPRINT",  "EPS",  and  "NA"  in  214,  E6.1,I4  format. 
The  "ITERM"  variable  denotes  the  maximum  nunber  of  cycles  to  be 
performed  in  direct  search.  "IPRINT"  is  a  print  control  para¬ 
meter  which  determines  if  printing  of  intermediate  optimization 
results  is  included  in  the  output,  The  "EPS"  parameter  is  the 
exit  tolerance  for  convergence  of  the  optimized  parameters  and 
"NA"  denotes  the  number  of  mesh  points  for  the  outer  region 
numerical  solution  to  the  similarity  equation.  A  typical  value 
of  "NA"  is  350  to  ensure  good  accuracy  and  should  be  increased 
in  applications  where  there  is  an  Intense  velocity  variation 
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in  the  outer  layer.  After  reading  in  all  of  these  variables, 
the  program  then  calls  "SUBROUTINE  PROFIT"  which  reads  in  exper¬ 
imental  data  pertinent  to  the  run. 

"SUBROUTINE  PROFIT"  begins  by  reading  in  a  four  digit 
identification  number  for  the  data  run  "ID".  The  choice  of 
an  identifier  is  up  to  the  user^s  discretion  and  will  be  suc¬ 
cessively  incremented  by  one  for  each  subsequent  data  station 
that  is  examined  in  the  data  set  under  consideration.  The  "ID" 
identifier  is  followed  by  "lUINT"  which  denotes  whether  the 
input  experimental  data  is  in  English  (ID*0)  or  metric  units 
(ID=1),  The  selection  of  metric  or  English  units  is  important 
to  the  input  since  this  determines  a  conversion  factor  which 
is  used  in  the  code.  A  user  supplied  title  for  the  data  set 
is  then  read  which  is  printe*  in  the  output.  The  number  of 
data  stations  "NSTA"  is  read  followed  by  input  parameters  for 
plotting.  These  plot  paraiieters  are  read  as  "NCYC,  NPLTPP, 

XLC,  XL,  NL,  INCPL"  and  are  accompanied  by  self-explanatory 
comments  describing  their  function  in  the  code.  Their  purpose 
is  to  allow  maximum  flexibility  in  obtaining  plotted  profiles 
to  meet  the  requirements  of  the  user.  The  only  plotting  comnent 
worth  mentioning  here  is  that  the  "NPLTPP"  parameter  denotes 
the  nunber  of  curves  to  be  plotted  on  one  plot  using  the  shifted 
Y-axis  method.  If  there  are  more  data  stations  than  the  NPLTPP 
parameter  CNSTA>NPLTPP),  then  another  plot  will  be  drawn  to 
plot  the  remaining  curves. 


The  data  that  follows  is  read  for  each  data  station  and 
the  sequence  described  in  this  paragraph  is  repetitive  for 
subsequent  data  stations.  First,  the  kinematic  viscosity  "XNU" 
is  read  in  E10.3  format  followed  by  "X,  UE,  UTAUE,  DELTA,  DELSTR, 
THETA,  DUEDX,  BETE”  in  7F10.0,  F10.2  format;  these  variables 
represent  the  x-location  on  the  wall,  local  mainstream  velocity, 
experimental  value  of  u^,  experimental  value  of  the  boundary 
layer  thickness  6,  experimental  value  of  displacement  thick¬ 
ness  experimental  value  of  momentum  thickness  e,  local 

value  of  the  velocity  gradient  dU^/dx,  and  the  value  of  tne 
Clauser  pressure  gradient  parameter  respectively.  An  optional 
read  statement  follows  which  reads  the  value  of  mainstream 
turbulence  level  "TUT”  for  examination  of  mainstream  turbulence 
effects.  The  number  of  experimental  velocity  data  points  "NDP" 
is  read  next  followed  by  the  starting  values  of  the  model  para¬ 
meters  K,  S,  and  <  denoted  by  xi^(l),  x^(2),  and  xi^(3)  respec¬ 
tively.  These  starting  values  are  used  to  initiate  optimiza¬ 
tion  in  the  direct  search  subroutine.  The  experimental  velocity 
profile  data  points  are  then  read  in  FI 0.0  format  starting  with 
the  x-coordinate  values  which  signify  the  distance  from  the 
wall  "TO”.  These  are  followed  by  the  corresponding  y-coordi nates 
values  which  denote  the  nondimensionalized  velocities  ”UD” 
represented  by  the  ratio  of  velocity  over  the  local  mainstream 
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velocity  u/Ug.  "SUBROUTINE  PROFIT"  then  calls  "SUBROUTINE  DRSRCH" 
to  initiate  the  optimization  of  the  profile  parameters. 

The  direct  search  subroutine  calls  all  other  subroutines 
to  evaluate  the  optimization  function  F(x).  After  the  optimiza¬ 
tion  has  been  completed,  information  pertinent  to  the  model  is 
printed  out  and  a  plot  of  the  analytical  profile  and  experi¬ 
mental  data  points  is  made.  A  listing  of  the  output  for  one 
data  station  is  presented  on  the  subsequent  pages.  An  output 
sunmary  similar  to  that  in  Table  5.1  for  all  data  stations  is 
printed  at  the  end  of  the  output  listings  for  all  data  stations. 


UT  MAINSTREAM  TUReULENCE  TeST  DATA 


1001 


X  a  40.300  FT  UEa  98.76  FT/SEC 


EXPERIMENTAL  VALUES 


OUEOXs 

3.000  1/SEC 

UTAU= 

4.026  FT/SEC 

DELTA* 

.6600  IN 

DELTA** 

.0824  IN 

THETA* 

.0673  IN 

EXPERIMENTAL 

VELOCITY  PROFILE 

61 

PO INTS 

Y  :iN) 

U/UE 

0. 0000 

0.0000 

.0063 

.3870 

.0  066 

.4110 

.0076 

.  4450 

.0086 

.4760 

.0066 

.4600 

.0111 

.5230 

.0123 

.6340 

.0133 

.6460 

.0166 

.6610 

.0176 

.5730 

.0167 

.6830 

.0211 

.6600 

.0223 

.6640 

.0247 

.6010 

.0267 

.6090 

.0287 

.6160 

.0302 

.6200 

.0  368 

.6360 

.0438 

.6600 

.06  06 

.6660 

.0566 

.6760 

.0633 

.6840 

.0707 

.6660 

.0766 

.7030 

.0830 

.7140 

.0607 

.7230 

.0663 

.7290 

.1034 

.7300 

.1103 

.7450 

.1163 

.7510 

.1237 

.7590 

.1308 

.7660 

.1473 

.7830 

.1661 

.8010 

.1824 

.8160 

.2007 

.8300 

.2178 

.8460 

.2363 

.8600 

.2628 

.8720 

.27  06 

.8860 

.2877 

.8600 

.3069 

.6110 

.3406 

.9330 

.3766 

.6620 

.4106 

.6680 

.4466 

.6020 

.4806 

.6610 
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8r'-ni<aaiM<yi  ipjii  «> 


A 


.  ?iJ3 

.  9960 

.5505 

.9980 

.5053 

1.0000 

.6205 

1.0000 

.6553 

1.0000 

.6906 

1.0000 

.7255 

1.0000 

.7604 

1.0000 

.7953 

.999  0 

.8309 

1.0000 

.  8653 

1.0000 

.9003 

.9990 

.9355 

.9990 

CYCLE 

c 

Xltl 

E 

0 

.1Z871ZE-03 

.168830C-ai 

F 

1 

.103S42F-03 

.1780a0E-01 

z 

.101538E-03 

.iasooac-01 

u 

3 

.  il73Z8E-fl3 

.198008E-01 

I 

k 

.100118e-03 

.ia4343E-01 

I 

5 

.100113E-e3 

.184ZtZE-01 

I 

If* 

.10011ZE>03 

.i84iiac-ai 

I 

7 

.lOailZE-03 

.184131E-01 

0*«IJ 

. lOOOOaE-02 
.IBDOaiE-OZ 
.loasooE-oz 
.  tOOOOOE-OZ 
. IIOOOOE-OZ 
.lOoatiE'OZ 
.laoossE-oz 
. taaosa e-oz 


r  ITE94TI0HS  IN  OBSHCH,  flX)»  .lOeBSEE-ai  NIH  *  1 
CAPK  >  .laitiaiE-ai  s  >  .iiez54E*oz  kapp*  >  .MraeaE^ao 


CIN>  .573116E*01 


TaPLUS>  .Z761SSE*OZ 


SETAC  X  a. 


NATCH  CONOITIOH  CRITEBIA  FOB  VEUOClTf  PROFILE 
ETAH  a  .0Z06  N«  >  390  ETAO  ■  .9Z0  8 


length  scales 

OELIN  -  .00004  OCLOUr  «  .16300 


PARAHCIEB  FIT 
UTAU  «  .4198B9E*a 
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» 

» 

EXPfRIweNTAL 

PROFILE  • 

• 

ANALYTICAL 

PROFILE  • 

U,  DEVIATION 

Y*  * 

UO/UE 

00*  * 

UA/UE 

UA*  » 

00*  -  UA» 

n 

a. 00  • 

O.OOOQ 

0.00  • 

O.OODD 

0.80  • 

Q.OOOOOr 

10.98  • 

,3870 

9.19  • 

.3458 

8.21  • 

.  977  821 

13,47  • 

.4110 

9.76  • 

.3926 

9.32  • 

.436  896 

15.<!4  • 

.445  0 

10.57  • 

.4249 

10.09  * 

.476  712 

17.61  • 

.475  0 

11.28  » 

.4521 

10.74  • 

.543  342 

19.68  • 

.4900 

11.64  • 

.4750 

11.28  * 

.356673 

a2.99  • 

.5230 

12.42  • 

.5043 

11.98  • 

.444  207 

25.48  • 

.5340 

12.68  • 

.5217 

12.39  * 

.292  495 

27.55  • 

.5460 

12.97  • 

.5339 

12.58  • 

.288  464 

32.11  * 

.6610 

13.32  • 

.5552 

13.18  • 

.138  810 

36.46  • 

.573  0 

13.61  • 

.5707 

13.55  • 

.055549 

40.81  • 

.5830 

13.84  • 

.5832 

13.85  • 

-.004265 

43.71  • 

.6900 

14.01  • 

.5904 

14.02  * 

-.009292 

46,20  • 

.5940 

14.11  ♦ 

.5960 

14.15  • 

-.0  4  8  205 

51.17  • 

.6010 

14.27  • 

.6062 

14,39  • 

-.122776 

55.31  • 

.6090 

14.46  • 

.6137 

14.57  * 

-.111887 

59.45  • 

.6150 

14.60  • 

.6206 

14.74  • 

-.133214 

€2.56  • 

.6200 

14,72  • 

.6254 

14.85  • 

-.128945 

76.23  • 

.6350 

15.08  • 

.6439 

15.29  • 

-.210444 

90.73  * 

.6500 

15.44  • 

.6598 

15.67  * 

-.233038 

105,24  • 

.665  0 

15,79  • 

.6732 

15.99  • 

-.194993 

117.25  • 

.675  0 

16.03  * 

.6829 

16.22  • 

-.187  475 

131.13  • 

.6840 

16.24  • 

.6929 

16.45  • 

-.211514 

146.46  • 

.695  0 

16.50  • 

.7029 

16.69  » 

-.188  439 

158.48  • 

.7030 

16.69  • 

.7103 

16.87  • 

-.172  684 

173.60  • 

,7140 

16,96  • 

.7191 

17.08  * 

-.119974 

187,89  • 

.7230 

17.17  • 

.7270 

17.26  » 

-.095266 

199.49  • 

.7290 

17,31  • 

.7333 

17.41  • 

-.101808 

214.20  • 

.738  0 

17,53  • 

.7411 

17.60  • 

-.072509 

228.49  • 

.7450 

17,69  • 

.7  484 

17.77  • 

-.061438 

240.92  • 

.7510 

17.83  • 

.7547 

17,92  • 

-.046436 

256.25  • 

.7590 

18.02  • 

.7624 

18.10  » 

-.079563 

270.96  • 

.766  0 

18.19  • 

.7695 

18.27  * 

-.083  912 

305.14  • 

.7830 

18,59  • 

.7858 

18.66  » 

-.065391 

342.02  • 

.8010 

19.02  • 

.8025 

19.06  • 

-.036515 

377.85  • 

.8160 

19.38  • 

.8181 

19.43  • 

-.0  5  0  706 

415.76  • 

.8300 

19,71  • 

.8339 

19.80  * 

-.091610 

451.19  • 

.8460 

20.09  • 

.8  478 

20.13  • 

-.042955 

487.44  • 

.8600 

20,42  * 

.  8613 

20.46  • 

-.031608 

523.69  • 

.8720 

20.71  * 

.8741 

20.76  » 

-.049303 

561.19  * 

.8060 

21.04  * 

.8864 

21.06  * 

-.0  1  0  40  8 

595.99  • 

.8980 

21.32  * 

.8972 

21.30  • 

.019704 

633.69  * 

.9110 

21.63  • 

.9080 

21.56  • 

.071655 

705.37  • 

.933  0 

22.16  • 

.9263 

22.00  • 

.159205 

778.08  • 

,952  0 

22.61  * 

.9420 

22.37  • 

.237851 

851.21  • 

.968  0 

22.99  • 

.9551 

22.68  • 

.307  042 

922.89  • 

.982  0 

23.32  • 

.9655 

22.93  • 

.390  605 

995.60  • 
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